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ABSTRACT 

We  report  parallax  measurements  for  70  ultracool  dwarfs  (UCDs)  including  11  late-M,  32  L,  and  27  T  dwarfs.  In 
this  sample,  14  M  and  L  dwarfs  exhibit  low  surface  gravity  features,  6  are  close  binary  systems,  and  2  are  metal-poor 
subdwarfs.  We  combined  our  new  measurements  with  114  previously  published  UCD  parallaxes  and  optical-mid- 
IR  photometry  to  examine  trends  in  spectral- type/ absolute  magnitude,  and  color-color  diagrams.  We  report  new 
polynomial  relations  between  spectral  type  and  MJHK.  Including  resolved  L/T  transition  binaries  in  the  relations, 
we  find  no  reason  to  differentiate  between  a  “bright”  (unresolved  binary)  and  a  “faint”  (single  source)  sample  across 
the  L/T  boundary.  Isolating  early  T  dwarfs,  we  find  that  the  brightening  of  T0-T4  sources  is  prominent  in  Mj  where 
there  is  a  [1.2-1. 4]  mag  difference.  A  similar  yet  dampened  brightening  of  [0.3-0. 5]  mag  happens  at  MH  and  a 
plateau  or  dimming  of  [—0.2  to  —0.3]  mag  is  seen  in  MK.  Comparison  with  evolutionary  models  that  vary  gravity, 
metallicity,  and  cloud  thickness  verifies  that  for  L  into  T  dwarfs,  decreasing  cloud  thickness  reproduces  brown 
dwarf  near-IR  color-magnitude  diagrams.  However  we  find  that  a  near  constant  temperature  of  1200  ±100  K  along 
a  narrow  spectral  subtype  of  T0-T4  is  required  to  account  for  the  brightening  and  color-magnitude  diagram  of  the 
L-dwarf/T-dwarf  transition.  There  is  a  significant  population  of  both  L  and  T  dwarfs  which  are  red  or  potentially 
“ultra-cloudy”  compared  to  the  models,  many  of  which  are  known  to  be  young  indicating  a  correlation  between 
enhanced  photospheric  dust  and  youth.  For  the  low  surface  gravity  or  young  companion  L  dwarfs  we  find  that  8  out 
of  10  are  at  least  [0.2-1. 0]  mag  underluminous  in  MJH  and/or  MK  compared  to  equivalent  spectral  type  objects. 
We  speculate  that  this  is  a  consequence  of  increased  dust  opacity  and  conclude  that  low  surface  gravity  L  dwarfs 
require  a  completely  new  spectral- type/ absolute  magnitude  polynomial  for  analysis. 
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1.  INTRODUCTION 

For  any  new  class  of  astronomical  objects,  distances  are  cru¬ 
cial  for  investigating  basic  physical  properties.  Brown  dwarfs, 
low-mass  objects  that  lack  sustained  stable  hydrogen  burning  in 
their  cores,  are  a  recent  addition  to  the  plethora  of  objects  studied 
in  astronomy.  They  were  first  predicted  by  Kumar  (1962)  and 
Hayashi  &  Nakano  (1963)  but  not  observationally  confirmed 
until  the  late  1990s  (Nakajima  et  al.  1995;  Rebolo  et  al.  1995). 
They  have  masses  between  ~0.072  M©  and  0.012  M©,  strad¬ 
dling  the  boundary  between  the  lowest  mass  stars  and  the  highest 
mass  exoplanets  (Saumon  et  al.  1996;  Chabrier  &  Baraffe  1997). 
In  early  2000,  the  standard  stellar  spectral  classification  scheme 
was  extended  beyond  M  dwarfs  to  include  “L”  dwarfs,  objects 
with  temperatures  ranging  between  1300  and  2000  K  and  “T-Y” 
dwarfs,  objects  cooler  than  1300  K  (see  Kirkpatrick  2005  and 
references  therein;  Cushing  et  al.  2011).  The  L-spectral  class 
encompasses  both  low-mass  stars  and  brown  dwarfs,  depending 
on  the  mass/age  of  the  ultracool  dwarf  (UCD;  see  Burrows  et  al. 
2001,  and  references  therein). 

9  Visiting  astronomer,  Cerro  Tololo  Inter- American  Observatory,  National 
Optical  Astronomy  Observatory,  which  are  operated  by  the  Association  of 
Universities  for  Research  in  Astronomy,  under  contract  with  the  National 
Science  Foundation. 
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Distances  provide  a  direct  means  for  calculating  the  luminos¬ 
ity  and  (if  there  is  a  reliable  radius  estimate)  the  effective  tem¬ 
perature  of  a  star  or  brown  dwarf.  Using  luminosity,  color,  and 
a  probe  of  effective  temperature,  brown  dwarf  color-magnitude 
diagrams  can  be  populated  and  used  to  investigate  physical  and 
chemical  parameters  such  as  gravity,  metallicity,  or  dust  proper¬ 
ties.  Moreover,  parallaxes  to  a  significant  number  of  objects  are 
required  to  create  relations  with  spectral  type  that  can  be  used 
to  estimate  distances  to  the  majority  of  brown  dwarfs  lacking 
astrometric  measurements. 

There  are  now  nearly  1000  spectroscopically  confirmed  field 
L  and  T  dwarfs11  that  define  observational  near-IR  color  trends. 
Spectroscopic  standards  have  been  designated  at  each  subtype 
that  exhibit  characteristic  features  of  the  changing  brown  dwarf 
spectral  energy  distribution.  The  physical  parameter  that  drives 
the  major  changes  in  the  observable  photometric  and  spectro¬ 
scopic  features  of  the  brown  dwarf  population  is  a  decreasing 
effective  temperature  (Teff).  However,  with  an  ever-growing  list 
of  brown  dwarfs  observed  in  the  field,  a  number  of  outliers 
have  emerged  that  exhibit  how  secondary  parameters  such  as 
age,  metallicity,  and  cloud  variability  can  change  observable 
properties.  There  are  a  handful  of  brown  dwarfs  that  have  halo 


11  According  to  the  dwarfarchives  Web  site  maintained  at 
http :  //  dwarfarchives .  org . 
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kinematics,  exhibit  blue  near-IR  colors,  and  have  enhanced 
metal  hydride  bands  along  with  weakened  metal  oxide  absorp¬ 
tion  bands  indicating  that  they  are  old  and  metal-poor  subdwarfs 
(Burgasser  et  al.  2003a,  2007;  Cushing  et  al.  2009;  Kirkpatrick 
et  al.  2010).  A  number  of  objects  have  red  near-IR  colors,  weak 
alkali  lines,  enhanced  metal  oxide  absorption  bands  in  the  opti¬ 
cal,  and  appear  to  be  juvenile-aged  members  of  nearby  moving 
groups  such  as  AB  Doradus,  /3  Pictoris,  etc  (Kirkpatrick  et  al. 
2006, 2010;  Cruz  et  al.  2009;  Allers  et  al.  2010;  Rice  et  al.  2010). 
There  are  also  objects  that  do  not  exhibit  the  extreme  spectral 
features  of  subdwarfs  or  low-gravity  dwarfs,  but  are  neverthe¬ 
less  near-IR  photometric  outliers  whose  photometric  properties 
might  be  attributable  to  dust,  subtle  variations  in  age  or  metal- 
licity,  or  photometric  variability  (Knapp  et  al.  2004;  Burgasser 
et  al.  2008b;  Faherty  et  al.  2009;  Schmidt  et  al.  2010;  Kirkpatrick 
et  al.  2010;  Looper  et  al.  2008b,  Radigan  et  al.  2012). 

With  the  early  parallax  programs  of  Dahn  et  al.  (2002), 
Tinney  et  al.  (2003),  and  Vrba  et  al.  (2004)  as  well  as  subsequent 
astrometric  programs,  the  spectral-type/absolute  magnitude 
relations  and  color-magnitude  diagrams  for  brown  dwarfs  were 
first  investigated  (Thorstensen  &  Kirkpatrick  2003;  Costa  et  al. 
2006;  Burgasser  et  al.  2008c;  Lucas  et  al.  2010;  Artigau  et  al. 
2010;  Schilbach  et  al.  2009;  Marocco  et  al.  2010).  One  of 
the  remarkable  features  of  UCD  color-magnitude  diagrams 
is  the  significant  scatter  in  luminosity  found  among  objects 
with  similar  spectral  types.  Increasingly  complex  atmospheric 
and  evolutionary  models  have  explained  this  as  the  result  of 
variations  in  secondary  parameters  such  as  gravity,  metallicity, 
sedimentation  efficiency,  and/or  binarity  (e.g.,  Tsuji  et  al.  1996; 
Tsuji  &  Nakajima  2003;  Burrows  et  al.  2006;  Helling  et  al.  2008; 
Saumon  &  Marley  2008).  The  models  disagree  as  to  which 
parameter  has  the  largest  effect  on  the  emergent  spectra  and 
color  trends,  and  small  numbers  of  objects  with  independently 
measured  secondary  parameters  have  hindered  progress  (Knapp 
et  al.  2004;  Patten  et  al.  2006;  Leggett  et  al.  2010).  Only 
by  increasing  the  number  of  well-characterized  UCDs  with 
distance  measurements  can  we  hope  to  understand  the  source 
and  extent  of  the  variation  in  the  color-magnitude  and  H-R 
diagrams. 

Another  major  anomaly  of  brown  dwarf  color-magnitude 
diagrams  is  an  intriguing  brightening  (or  bump)  in  the  J  band 
(up  to  1.5  mag;  Vrba  et  al.  2004;  Tinney  et  al.  2003),  and 
to  some  extent  H  and  K ,  as  objects  transition  between  the 
warmer  L  dwarf  and  cooler  T  dwarf  spectral  classes.  In  the 
past  decade,  several  L/T  transition  objects  have  been  confirmed 
as  flux  reversal  binaries  with  resolved  components  straddling  the 
bump.  These  objects  confirm  that  the  brightening  is  an  intrinsic 
feature  of  brown  dwarf  evolution  (e.g.,  Looper  et  al.  2008a;  Gizis 
et  al.  2003;  Liu  et  al.  2006;  Burgasser  et  al.  2006b;  Burgasser 
2007);  however,  the  full  extent  and  physical  explanation  remain 
mysterious.  One  possible  explanation  is  that  opacity  from 
condensate  clouds  is  especially  influential  on  the  1  fim  region  of 
L  and  early-type  T  dwarfs,  so  changes  in  cloud  properties  can 
explain  the  unusual  brightening  (e.g.,  Burgasser  et  al.  2002b; 
Knapp  et  al.  2004).  The  small  numbers  of  L/T  transition 
objects  with  parallax  measurements  has  hindered  progress  in 
understanding  the  degree  and  variation  in  brightening  exhibited 
at  this  interesting  phase  of  brown  dwarf  evolution. 

In  late  2006,  we  initiated  the  Brown  Dwarf  Kinematics  Project 
(BDKP)  in  order  to  address  persistent  questions  of  brown  dwarf 
evolution  and  atmospheric  properties  using  astrometric  mea¬ 
surements  of  proper  motion,  parallax,  and  radial  velocity.  Proper 
motion  analysis  of  the  population  was  reported  in  Faherty  et  al. 


(2009,  2010,  2011).  In  this  work,  we  report  parallaxes  for 
70  UCDs.  Section  2  describes  the  target  list  as  well  as  the 
data  acquisition  and  reduction.  Section  3  describes  the  paral¬ 
lax  pipeline  used  to  determine  distances.  Section  4  uses  all 
parallax  measurements  reported  in  this  work  in  combination 
with  published  values  and  photometric  information  obtained 
from  various  catalogs  to  investigate  spectral  type/absolute  mag¬ 
nitude  relations  and  color-magnitude  diagrams  for  the  brown 
dwarf  population.  In  Section  5  an  updated  brown  dwarf  near- 
IR  color-magnitude  diagram  is  examined  using  evolutionary 
models.  Sections  6  and  7  discuss  the  absolute  magnitude  versus 
spectral  type  relation  for  low  surface  gravity  and  subdwarfs, 
respectively.  Section  8  reviews  the  kinematics  for  an  ensem¬ 
ble  of  all  known  brown  dwarfs  with  parallax  measurements  and 
Section  9  investigates  both  known  and  suspected  binaries  within 
the  sample.  Conclusions  are  reported  in  Section  10. 

2.  OBSERVATIONS 
2.7.  Target  List 

We  compiled  the  parallax  target  list  from  the  BDKP  astro- 
metric  sample  reported  in  Faherty  et  al.  (2009).  Instrumental 
limitations  precluded  measuring  parallaxes  to  the  faintest,  most 
distant  L  and  T  dwarfs  so  we  focused  primarily  on  objects 
within  20  pc.  However,  we  were  also  interested  in  subsets  of  the 
population  which  included  low  surface  gravity  dwarfs  (poten¬ 
tially  young  sources)  and  subdwarfs  (potentially  old  sources). 
For  these  scientifically  interesting  subsets  we  relaxed  our  astro- 
metric  constraint  to  include  sources  whose  predicted  spectro¬ 
scopic  parallax  was  up  to  50  pc.  Our  full  target  list  consisted  of 
70  dwarfs,  including  1 1  M,  32  L,  and  27  T  dwarfs  (see  Table  1). 
In  this  sample  there  were  14  low  surface  gravity  dwarfs,  2  subd¬ 
warfs,  6  binary  systems,  and  9  calibrators  with  previous  parallax 
measurements. 

2.2.  Data  Collection  and  Reduction 
2.2.1.  ANDICAM 

We  obtained  parallax  frames  with  the  ANDICAM  (A  Novel 
Double-Imaging  CAMera-  DePoy  2003)  imager  between  2006 
November  and  2010  March  (~500  hr  of  observations).  All  data 
were  acquired  through  queue  observing  with  telescope  time  al¬ 
located  through  the  SMARTS  (Small  and  Moderate  Aperture 
Research  Telescope  System)  consortium.  ANDICAM  is  a  dual 
channel  near-IR  and  CCD  imager  located  on  the  1 .3  m  telescope 
at  Cerro  Tololo  Inter  American  Observatory  (CTIO).  The  optical 
detector  is  a  1024  x  1024  CCD  and  the  near-IR  channel  uses 
a  Rockwell  1024  x  1024  HgCdTe  HAWAII  array.  The  near-IR 
field  of  view  is  2.4  arcmin  with  an  unbinned  plate  scale  of 
0!  137  pixel-1.  The  optical  CCD  field  of  view  is  larger, 
6.2  arcmin,  with  a  plate  scale  of  07369  pixel-1.  The  optical 
and  near-IR  channels  operate  independently  with  a  dichroic  fil¬ 
ter  directing  light  to  the  two  independent  cameras.  Therefore, 
we  were  able  to  take  a  set  of  near-IR  images  while  integrating 
in  the  optical. 

To  ensure  the  same  reference  stars  for  each  parallax  frame, 
we  required  the  target  star  to  always  be  placed  in  the  same  X, 
Y  position  on  the  detector.  We  also  required  all  observations 
to  be  made  within  ±30  minutes  of  meridian  crossing  to 
minimize  the  corrections  for  differential  color  refraction  (DCR, 
see  Section  3.1  below).  Typical  seeing  was  1"  and  useable 
conditions  for  our  parallax  program  were  up  to  2!' . 
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Table  1 

Target  List 


Name 

SpT 

SpT 

Nights 

Framesa 

Ref  Starsa 

At 

J 

Instrument*5 

Notec 

Flagd 

Ref 

OPT 

IR 

(yr) 

(MKO) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

2MASSJ00325584— 4405058 

LO.Oy 

- 

12 

49 

71 

2.24 

14.68  ±  0.04 

A 

LG 

23,22 

2M  ASS  J00345 1 57+0523050 

- 

T6.5 

6 

27 

38 

2.91 

15.11  ±0.03 

I 

F 

24,  19 

2MASSIJ0103320+193536 

L6.O0 

- 

6 

30 

54 

3.08 

16.16  ±0.08 

I 

LG 

F 

15 

2MASSJ02212859— 6831400 

M8.O0 

- 

15 

60 

78 

2.04 

13.89  ±  0.03 

A 

LG 

17 

2MASSJ03 185403—3421292 

L7.0 

- 

8 

34 

46 

3.08 

15.44  ±  0.05 

I 

F 

16 

SDSS  J032553 . 17+042540. 1 

- 

T5.5 

6 

33 

21 

3.08 

15.90  ±  0.03 

I 

F 

2 

2MASSJ0334121 8—4953322 

M9.0 

- 

16 

54 

73 

2.96 

11.32  ±0.02 

A 

77 

2MASSJ04221413+1530525: 

M6.0y 

- 

11 

39 

38 

2.97 

12.67  ±  0.02 

A 

LG 

17 

2MASSJ04390101 -2353083 

L6.5 

- 

14 

49 

58 

3.00 

14.31  ±0.03 

A 

3 

2MASSJ04455387— 3048204 

L2.0 

- 

25 

84 

71 

3.05 

13.31  ±0.03 

A 

3 

2M  AS  S  J05 0 1 2406 — 00 1 045 2 

L4.0y 

- 

11 

46 

89 

2.92 

14.86  ±  0.04 

A 

LG 

17,  22, 

2MASSJ05160945— 0445499 

- 

T5.5 

11 

61 

34 

3.90 

15.95  ±  0.08 

I 

26,19 

2MASSJ05 184616—2756457 

LO.Oy 

- 

10 

42 

42 

3.14 

15.16  ±0.04 

I 

LG 

N 

13 

2MASSJ05 185995—2828372 

L7.5 

T1.0 

11 

55 

66 

3.90 

15.87  ±0.10 

I 

B 

27,  16,19 

2MASSJ05361998- 1920396 

LI. 00 

- 

9 

38 

59 

3.14 

15.65  ±  0.08 

I 

LG 

N 

13 

LHS1777 

M5.0 

- 

21 

79 

32 

2.84 

10.15  ±0.02 

A 

Cal 

78 

2MASSJ06085283— 2753583 

M8.5y 

- 

22 

80 

92 

3.04 

13.53  ±  0.03 

A 

LG 

3,79 

2M  AS  S  J06 1 64006 — 6407 1 94 

sdL5.0 

sdL5.0 

8 

28 

78 

3.06 

16.40  ±0.11 

I 

SD 

F,N 

30 

2MASSJ06244595— 4521548 

L5.0 

- 

17 

63 

86 

2.27 

14.36  ±  0.03 

A 

17 

DENIS-PJ065248. 5— 574137 

M8.O0 

- 

12 

52 

95 

2.83 

13.56  ±  0.03 

A 

LG 

17 

2MASSJ07123786— 6155528 

LI. 00 

- 

8 

36 

88 

3.14 

15.20  ±  0.06 

I 

LG 

F,N 

22 

2MASSJ07290002— 3954043 

- 

T8.0 

11 

39 

74 

3.70 

15.88  ±0.08 

I 

25 

SDSSJ074201. 41+205520.5 

- 

T5.0 

7 

19 

75 

2.08 

15.60  ±  0.03 

I 

F 

9,19 

2MASSIJ0746425+200032 

L0.5 

L1.0 

13 

52 

91 

2.69 

11.64  ±0.03 

A 

B,Cal 

10,  15,9 

SDSSJ083048. 80+01283 1 . 1 

- 

T4.5 

9 

51 

46 

3.90 

15.99  ±  0.03 

I 

9,19 

2MASSJ08472872- 1532372 

L2.0 

- 

24 

85 

69 

2.89 

13.42  ±  0.03 

A 

3 

2MASSIJ0859254- 194926 

L7.0 

- 

9 

38 

32 

3.71 

15.42  ±  0.05 

I 

3 

2MASSJ09393548— 2448279 

- 

T8.0 

8 

33 

63 

3.81 

15.61  ±  0.09 

I 

B,Cal 

F 

31,19 

2MASSJ09490860- 1545485 

- 

T2.0 

11 

61 

25 

3.90 

16.09  ±0.12 

I 

31,19 

WT248 

M3.0 

- 

15 

63 

99 

2.06 

10.51  ±0.02 

A 

Cal 

80 

2MASSJ10073369— 4555147 

- 

T5.0 

8 

30 

97 

3.70 

15.65  ±  0.07 

I 

F 

25 

2MASSIJ1010148— 040649 

L7.0 

- 

10 

45 

55 

3.70 

15.39  ±  0.06 

I 

3 

2MASSJ1 0220489+0200477 

M9.O0 

- 

6 

34 

57 

3.06 

14.03  ±  0.03 

I 

LG 

F,N 

17 

SDSSJ103026.78+021 306.4 

- 

L9.5 

8 

41 

52 

3.70 

17.10  ±0.05 

I 

F,N 

9 

2MASSWJ1036530— 344138 

L6.0 

- 

6 

29 

69 

2.82 

15.51  ±0.05 

I 

F 

32 

SDSSJ104335. 08+121 3 14.1 

- 

L7.0 

7 

30 

63 

3.71 

15.82  ±  0.03 

I 

F 

2 

SDSSJ1 04409.43+042937.6 

- 

L7.0 

8 

40 

27 

3.90 

15.84  ±  0.03 

I 

F 

9 

2MASSJ10584787- 1548172 

L3.0 

L3.0 

14 

56 

28 

2.48 

14.12  ±0.05 

A 

Cal 

33,  14,9 

TWA28 

M8.5y 

- 

18 

67 

97 

2.89 

12.96  ±  0.02 

A 

LG 

81 

2MASSJ1 1145133-2618235 

- 

T7.5 

9 

42 

42 

2.08 

15.52  ±  0.05 

I 

31,19 

TWA26 

M9.0y 

- 

15 

57 

90 

2.22 

12.61  ±  0.03 

A 

LG 

17 

2MASSJ1 1553952-3727350 

L2.0 

- 

21 

91 

97 

2.53 

12.73  ±  0.02 

A 

32 

SDSSJ1 15553.86+055957.5 

- 

L7.5 

7 

31 

49 

3.70 

15.63  ±  0.03 

I 

F 

9 

SDSSJ1 15700.50+061 105.2 

- 

T1.5 

7 

37 

23 

3.90 

17.09  ±  0.05 

I 

F,N 

9,19 

SDSSJ120747. 17+024424. 8 

L8.0 

TO.O 

7 

32 

61 

1.88 

15.38  ±0.03 

I 

F,N 

6,19 

2MASSJ12095613- 1004008 

T3.5 

T3.0 

6 

30 

27 

3.70 

15.55  ±  0.03 

I 

B 

F,N 

24,  16,19 

2MASSJ12154432— 3420591 

- 

T4.5 

8 

45 

53 

3.90 

16.24  ±0.13 

I 

F,N 

25 

2MASSJ135955 10— 4034582 

L1.0 

- 

15 

53 

98 

1.80 

13.58  ±0.03 

A 

17 

2MASSJ14044941— 3159329 

TO.O 

T2.5 

9 

35 

87 

3.90 

15.51  ±0.06 

I 

B 

25,  35 

2M  AS  S  J 1 4442067 — 20 1 9222 

sdM9.0 

- 

15 

63 

87 

2.87 

12.51  ±0.02 

A 

SD,Cal 

72,  82 

SDSSpJ14460060+002452.0 

L6.0 

L5.0 

7 

36 

52 

2.08 

15.56  ±  0.05 

I 

Cal 

F,N 

34,  6,9 

SDSSJ15041 1.63+102718.4 

- 

T7.0 

9 

39 

66 

2.08 

16.49  ±  0.03 

I 

2 

SDSSJ151 1 14.66+060742.9 

- 

TO.O 

8 

46 

43 

2.08 

15.83  ±  0.03 

I 

B 

F 

2 

SDSSJ152103. 24+013 142.7 

- 

T2.0 

6 

38 

41 

2.08 

16.06  ±  0.03 

I 

F 

9,19 

2MASSIJ1526140+204341 

L7.0 

- 

6 

23 

73 

1.88 

15.48  ±  0.05 

I 

F 

15 

2MASSJ16150413+1340079 

- 

T6.0 

8 

33 

86 

2.08 

16.32  ±  0.09 

I 

F 

25 

SDSSJ163022.92+081 822.0 

- 

T5.5 

6 

25 

73 

1.88 

16.18  ±0.03 

I 

F,N 

2 

2MASSJ1645221 1  —  1319516 

L1.5 

- 

12 

47 

99 

1.16 

12.37  ±  0.03 

A 

32 

2MASSJ17545447+1649196 

- 

T5.5 

5 

11 

98 

1.88 

15.79  ±  0.07 

I 

F 

1 

2MASSJ1 8283572-4849046 

- 

T5.5 

9 

39 

93 

1.88 

15.17  ±0.06 

I 

24,19 

2MASSJ193601 87—5502322 

L5.0 

- 

13 

37 

97 

1.31 

14.40  ±  0.04 

A 

17 

SDSS  J204749 .6 1  -07 1 8 1 8 .3 

- 

TO.O 

6 

22 

79 

1.34 

16.70  ±  0.03 

I 

F 

9,19 

SDSS  J205235 .3 1  - 1 60929.8 

- 

T1.0 

6 

13 

71 

1.88 

16.04  ±  0.03 

I 

B 

F,N 

2 

2MASSJ20575409— 0252302 

L1.5 

LI. 5 

20 

75 

93 

2.20 

13.04  ±0.02 

A 

3,7 

3 


The  Astrophysical  Journal,  752:56  (22pp),  2012  June  10 


Table  1 

(Continued) 


Faherty  et  al. 


Name 

SpT 

SpT 

Nights 

Framesa 

Ref  Starsa 

At 

J 

Instrument15 

Notec 

Flagd 

Ref 

OPT 

IR 

(yr) 

(MKO) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

2MASSJ21321 145+1341584 

L6.0 

- 

5 

19 

75 

1.20 

15.68  ±  0.06 

I 

B 

F,N 

13 

2MASSJ21513839— 4853542 

- 

T4.0 

7 

29 

72 

1.52 

15.71  ±  0.07 

I 

F 

38,19 

2MASSJ22282889— 43 10262 

- 

T6.0 

7 

29 

52 

3.47 

15.64  ±  0.07 

I 

F 

26,19 

2MASS  J23224684-  3 133231 

L0.0£ 

- 

9 

34 

66 

1.39 

13.50  ±0.03 

A 

LG 

17 

2MASSIJ2356547— 155310 

- 

T5.5 

7 

27 

46 

3.47 

15.48  ±  0.03 

I 

Cal 

F 

36,19 

2MASSJ23594034— 7335055 

- 

T5.5 

9 

42 

55 

3.47 

16.10  ±0.10 

I 

1 

Notes. 

a  Number  of  frames  and  reference  stars  used  in  the  parallax  solution. 
b  I:  ISPI  and  A:  ANDICAM. 

c  LG  is  a  low  surface  gravity  dwarf,  Cal  is  a  calibrator  ultracool  dwarf,  SD  is  an  ultracool  subdwarf,  B  is  a  tight  binary  unresolved  in  2MASS. 
d  F  indicates  an  object  with  <9  parallax  frames;  A  indicates  an  object  with  Mjhk  uncertainties  >0.5  mag  therefore;  unless  otherwise  noted,  it  indicates  that  the  target 
was  not  used  in  the  analysis  throughout  the  paper. 

References.  (1)  This  paper;  (2)  Chiu  et  al.  (2006);  (3)  Cruz  et  al.  2003;  (4)  Gizis  et  al.  2000;  (5)  Golimowski  et  al.  2004;  (6)  Hawley  et  al.  2002;  (7)  Kendall  et  al. 
2004;  (8)  Kirkpatrick  et  al.  1997;  (9)  Knapp  et  al.  2004;  (10)  Reid  et  al.  2000;  (11)  Ruiz  et  al.  1997;  (12)  Wilson  et  al.  2001;  (13)  Cruz  et  al.  2007;  (14)  Kirkpatrick 
et  al.  1999;  (15)  Kirkpatrick  et  al.  2000;  (16)  Kirkpatrick  et  al.  2008;  (17)  Reid  et  al.  2008;  (18)  Forveille  et  al.  2004;  (19)  Burgasser  et  al.  2006a;  (20)  Burgasser 
et  al.  2008a;  (21)  Stumpf  et  al.  2008;  (22)  Cruz  et  al.  2009;  (23)  EROS  Collaboration  et  al.  1999;  (24)  Burgasser  et  al.  2003d;  (25)  Looper  et  al.  2007;  (26)  Burgasser 
et  al.  2004;  (27)  Cruz  et  al.  2004;  (28)  Burgasser  et  al.  2000b;  (29)  Burgasser  et  al.  2003b;  (30)  Cushing  et  al.  2009;  (31)  Tinney  et  al.  2005;  (32)  Burgasser  et  al. 
2002a;  (33)  Delfosse  et  al.  1997;  (34)  Geballe  et  al.  2002;  (35)  Looper  et  al.  2008a;  (36)  Burgasser  et  al.  2002a;  (37)  Lepine  et  al.  2003a;  (38)  Ellis  et  al.  2005;  (39) 
Martin  et  al.  1999;  (40)  Thorstensen  &  Kirkpatrick  2003;  (41)  Burgasser  et  al.  2004;  (42)  Fan  et  al.  2000;  (43)  Burgasser  et  al.  2003a;  (44)  Scholz  et  al.  2003;  (45) 
Leggett  et  al.  2000;  (46)  Luhman  et  al.  2007;  (47)  Lodieu  et  al.  2007;  (48)  Biller  et  al.  2006;  (49)  Strauss  et  al.  1999;  (50)  Burgasser  et  al.  1999;  (51)  Artigau  et  al. 
2010;  (52)  Tsvetanov  et  al.  2000;  (53)  Nakajima  et  al.  1995;  (54)  Mugrauer  et  al.  2006;  (55)  Pinfield  et  al.  2008;  (56)  Burgasser  et  al.  2000a;  (57)  Goldman  et  al. 
2010;  (58)  Burningham  et  al.  2009;  (59)  Warren  et  al.  2007;  (60)  Delorme  et  al.  2008;  (61)  Burningham  et  al.  2008;  (62)  Lucas  et  al.  2010;  (63)  Schneider  et  al.  2002; 
(64)  Burgasser  et  al.  2007;  (65)  Kirkpatrick  et  al.  2001;  (66)  Kasper  et  al.  2007;  (67)  Burgasser  et  al.  2010b;  (68)  Costa  et  al.  2006;  (69)  Dahn  et  al.  2002;  (70)  Henry 
et  al.  2006;  (71)  Perryman  et  al.  1997;  (72)  Schilbach  et  al.  2009. 


In  the  optical  we  observed  in  the  Ic  band  with  integration 
times  that  ranged  from  265  s  for  our  brightest  targets  to  610  s 
for  our  faintest  (typical  S/N  >  100  for  all  targets).  In  the  near- 
IR  we  observed  in  the  J  band  with  integration  times  that  ranged 
from  20  s  with  5  co-adds  for  our  brightest  targets  to  130  s  with 
8  co-adds  for  our  faintest.  We  acquired  5-7  near-IR  images  in  a 
10"  dither  pattern. 

The  optical  ANDICAM  data  were  processed  with  overscan 
subtraction  and  flat-fielding,  prior  to  distribution.  Initially  we 
intended  to  use  the  near-IR  data  for  parallaxes  and  the  CCD 
data  as  a  check  on  the  astrometric  quality.  However,  we  quickly 
realized  that  the  optical  images  were  far  superior  to  the  near- 
IR,  which  were  plagued  with  imaging  artifacts,  an  occasional 
elongated  point  spread  function  (PSF),  and  a  smaller  field  of 
view  (therefore  fewer  reference  stars).  As  a  result  we  report 
parallaxes  in  this  paper  based  only  on  the  optical  imaging. 

2.2.2.  Infrared  Side  Port  Imager  ( ISPI) 

We  collected  parallax  data  for  our  faintest  targets  with  the 
Infrared  Side  Port  Imager  (ISPI)  on  the  CTIO  4  m  Blanco 
telescope  (van  der  Bliek  et  al.  2004).  ISPI  is  a  2048x2048 
HgCdTe  HAWAII-2  array  with  an  ~8  arcmin  field  of  view 
and  nominal  plate  scale  of  O'/ 303  pixel-1.  Observations  were 
conducted  over  a  period  of  just  under  2  yr  (from  early  2008 
March  through  late  2010  April)  on  15  observing  runs.  After  the 
spring  of  2010,  ISPI  was  replaced  on  the  Blanco  4  m  telescope 
with  the  NEWFIRM  infrared  detector.  However,  it  was  put  back 
on  the  telescope  in  October  of  201 1  and  we  obtained  three  more 
epochs  over  six  nights  of  observing  in  2011  November,  2012 
January,  and  2012  February.  These  new  parallax  frames  are  also 
included  in  this  work. 

As  opposed  to  ANDICAM,  ISPI  data  were  collected  clas¬ 
sically.  All  observations  were  carried  out  in  the  J  band 


under  seeing  conditions  up  to  2"  FWHM  with  typical  condi¬ 
tions  between  Of! 8  and  1"  1 .  Most  of  the  parallax  observations 
were  made  when  the  target  was  within  ±30  minutes  of  the 
meridian  to  minimize  the  corrections  for  DCR  (see  Section  3.1 
below).  However,  due  to  observational  constraints  (weather,  in¬ 
strument  issues,  etc.)  some  targets  were  observed  within  ±1  hr 
of  meridian  crossing. 

In  order  to  minimize  the  effects  of  distortion  and  to  ensure  the 
same  reference  stars  in  each  frame,  we  placed  the  target  star  on 
the  same  X,  Y pixel  position  for  each  parallax  frame.  On  the  first 
observing  run  for  a  target,  the  frame  was  initially  offset  from 
the  center  of  the  chip  to  avoid  the  four-quadrant  seam  along  the 
detector.  This  initial  frame  was  used  in  all  subsequent  observing 
runs  as  a  reference  for  determining  where  to  place  the  parallax 
star. 

Integration  times  were  set  by  the  magnitude  of  the  target  and 
the  conditions  at  the  telescope.  They  ranged  from  30  to  60  s  with 
5-10  co-adds  and  5-10  images  in  a  10"  dither  pattern  (typical 
S/N  >  100  for  all  targets).  Depending  on  the  weather  and 
seeing  conditions,  the  typical  integration  time  per  target  was 
15-40  minutes. 

Dark  frames  and  lights  on/off  dome  flats  were  obtained  at 
the  start  of  each  evening.  Reduction  procedures  were  based  on 
the  prescriptions  put  together  by  the  ISPI  team12  utilizing  a 
combination  of  IRAF  routines.  J- band  flats  were  created  by 
median  combining  the  lights  on  and  lights  off  images  then 
subtracting  the  two.  Bad  pixel  masks  were  created  from  a 
dome  flat  image.  Individual  parallax  frames  were  flat-fielded 
and  corrected  for  bad  pixels  with  the  calibration  images.  All 
images  were  flipped  to  orient  north  up  and  east  to  the  left  using 
the  IRAF  routine  osiris  in  the  cirred  package.  Finally,  the  IRAF 


12  http://www.ctio.noao.edu/instruments/ir_instruments/ispi/ 
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routine  xdimsum  was  used  to  perform  sky  subtractions  and  mask 
holes  from  bright  stars.13 

3.  PARALLAX  PIPELINE 
3.1.  Source  Extraction 

Once  all  images  were  reduced,  we  used  the  Carnegie  As¬ 
trometric  Planet  Search  (hereafter  ATPa)  software  to  extract 
all  point  sources  and  solve  for  relative  parallaxes  and  proper 
motions  (Boss  et  al.  2009).  Images  were  not  co-added;  rather, 
sources  were  extracted  on  every  image  at  every  dither  position. 

The  processing  of  each  night  of  observation  (or  epoch)  was 
treated  separately  and  consisted  of  extracting  the  position  of 
all  point  sources.  This  required  the  specification  of  (1)  a  guess 
of  the  FWHM  of  each  set  of  images,  (2)  the  plate  scale  for 
the  instrument,  and  (3)  a  high  quality  image  to  be  used  as  the 
template  used  to  generate  a  preliminary  list  of  sources  to  be 
extracted  on  all  other  frames  in  a  given  epoch. 

The  precise  centroids  of  the  stars  were  measured  by  binning 
the  stellar  profile  in  the  X  and  Y  directions  using  a  box 
of  ~2"  around  the  pixel  with  maximum  flux.  Each  one¬ 
dimensional  profile  was  precisely  centroided  by  finding  the 
zero  of  the  profile  convolution  using  Tukey’s  biweight  function 
(Tukey  1977).  The  width  of  this  function  depends  on  a  scale 
parameter  to  compute  the  centroid  which  was  varied  and  the 
average  taken  (FWHM-0.5  pixel,  FWHM-0.25  pixel,  FWHM, 
FWHM+0.25  pixel,  FWHM+0.5  pixel)  to  mitigate  systematic 
effects  caused  by  poor  sampling  of  the  stellar  PSF.  We  tested 
a  variety  of  centroiding  methods  using  different  profiles  but 
decided  that  the  convolution  approach  described  above  provided 
the  best  centroid  accuracy  and  robustness  (see  Boss  et  al. 
2009  for  further  details).  Once  all  sources  were  extracted  on 
all  the  frames  in  a  given  night,  the  10  stars  showing  the  best 
formal  centroid  uncertainties  were  used  to  define  a  preliminary 
reference  frame  and  their  standard  deviation  assigned  as  the 
precision  per  image. 

ATPa  generates  a  text  file  (from  now  on  referred  to  as  a  plate 
file)  containing  X,  Y  positions  with  corresponding  uncertainties 
(in  pixels),  and  a  rough  flux  measurement  of  all  the  successfully 
extracted  stars.  Plate  files  were  used  in  the  next  step  of  the 
processing  to  calibrate  the  field  distortions  and  measure  the 
motion  of  each  star  as  a  function  of  time.  Objects  displaying 
positional  uncertainties  larger  than  5  times  the  median  precision 
were  automatically  removed.  These  were  usually  spurious 
sources  caused  by  unfiltered  detector  artifacts  or  cosmic  rays. 
The  typical  centroiding  uncertainty  for  parallax  targets  in  ISPI 
and  ANDICAM  was  ~0.01  pixels. 

DCR  corrections  are  typically  required  because  the  parallax 
star  and  reference  stars  have  very  different  colors.  As  a  result, 
their  positions  shift  relative  to  one  another  due  to  different 
amounts  of  atmospheric  refraction.  The  effect  is  wavelength, 
weather,  and  zenith-distance  dependent.  Stone  (1996,  2002) 
presents  a  theoretical  method  for  determining  DCR  effects.  That 
work  demonstrates  that  by  maintaining  small  zenith  distances, 
DCR  effects  in  I  and  longer  wavelengths  (such  as  J)  are  minimal, 
typically  <  1  mas.  Similar  results  were  found  using  the  empirical 
methodology  proposed  by  Monet  et  al.  (1992).  The  low-mass 
star  optical  parallax  program  of  Jao  et  al.  (2005)  and  the  brown 
dwarf  optical  parallax  program  of  Dahn  et  al.  (2002)  also 


13  We  note  that  during  observations  the  primary  mirror  would  occasionally 
vibrate  causing  the  PSF  to  appear  elongated.  When  this  occurred  we  would 
halt  and  restart  an  integration.  The  problem  was  sporadic  but  did  not  affect  any 
of  the  final  images. 


found  negligible  I- band  DCR  corrections  as  did  the  near-IR  T 
dwarf  parallax  program  of  Tinney  et  al.  (2003).  Therefore,  DCR 
corrections  are  not  applied  to  the  positions  in  our  pipeline.  To 
ensure  that  even  this  small  effect  was  minimized,  we  observed 
targets  (with  few  exceptions)  within  ±30  minutes  of  meridian 
crossing. 

3.2.  Astrometric  Solution  and  Parallax 

Using  the  Astrometric  Iterative  Solution  included  in  the  ATPa 
package,  the  extracted  X,  Y  pixel  positions  were  mapped  to 
a  local  tangent  plane  in  order  to  solve  for  the  astrometric 
parameters.  The  highest  quality  image  obtained  for  a  given 
parallax  target  was  designated  as  the  initial  catalog  to  which 
everything  was  matched  in  the  first  astrometric  iteration.  The 
solution  needed  to  be  initialized  with  the  right  ascension 
(R.A.)  and  declination  (decl.)  coordinates  of  the  brightest  star 
in  the  field;  however,  the  coordinates  were  only  required  to 
compute  the  projected  paralactic  motion  in  the  tangent  plane. 
Therefore  a  precision  of  a  few  arcseconds  on  this  bright  star 
was  more  than  sufficient  to  initialize  the  solution.  In  the  first 
iteration,  the  position  of  the  stars  in  the  initial  catalog  were 
still  approximate,  so  only  a  linear  transformation  was  applied  to 
each  plate  to  correct  for  telescope  pointing,  field  rotation,  and 
changes  in  the  plate  scale.  This  matching  provided  the  position 
of  each  star  as  a  function  of  time  in  the  reference  frame  defined 
by  the  initial  catalog.  The  apparent  trajectory  of  each  star  was 
then  fit  to  a  linearized  astrometric  model 

x(t)  =  x0  +  tia  (t  -  t0)  -  Upa(t) ,  (1) 

yit)  =  y0  +  [i8(t  -  t0)  -  n ps(t) ,  (2) 

where  xo  and  yo  are  positional  offsets  at  the  first  epoch  of 
observation  f0;  fia  and  ii8  are  the  proper  motion  in  R.A.  and 
deck,  respectively;  II  is  the  parallax;  and  pa  and  p$  are  the 
parallax  factors  in  R.A.  and  deck,  respectively.  At  this  point,  all 
the  quantities  are  given  in  milliarc seconds  (mas)  and  the  time  t 
is  measured  in  years.  The  parallax  factors  are  computed  using 
the  Earth  geocenter  as  obtained  from  the  DE405  Ephemeris.14 
This  linearized  model  is  based  on  the  prescriptions  laid  out  in 
the  HIPPARCOS  (Perryman  et  al.  1997)  and  Tycho  Catalogues15 
(Hog  et  al.  2000). 

After  the  astrometric  solution  of  the  field  was  obtained,  a 
subset  of  well-behaved  stars  (rms  <  5  mas  and  at  least  4-5 
observations)  was  chosen  and  the  second  iteration  begun  using 
those  as  reference  sources.  The  whole  process  (crossmatching, 
field  distortion  fit,  and  astrometric  solution)  was  repeated  a 
number  of  times.  After  the  first  iteration,  ATPa  allowed  fitting 
more  detailed  field  distortion  corrections  using  higher  order 
polynomials  in  X  and  Y.  We  tested  up  to  third-order  polynomials 
but  these  yielded  negligible  improvement  in  both  ANDICAM 
and  ISPI  images.  Consequently,  we  ran  all  targets  using  a 
second-order  polynomial. 

We  set  the  number  of  total  iterations  to  3-5.  The  convergence 
was  monitored  by  checking  the  average  rms  of  the  well-behaved 
reference  sources.  This  iterative  process  was  automated  by  ATPa 
and  required  little  supervision.  The  selection  of  reference  stars 
was  done  automatically,  but  ATPa  allowed  the  user  to  flag 
problematic  reference  stars.  As  a  result,  after  a  first  solution 
was  obtained,  we  checked  the  final  catalog  and  re-ran  the  whole 

14  http://ssd.jpl.nasa.gov/ 

15  http://www.rssd.esa.int/SA/HIPPARCOS/docs/voll_all.pdf 
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Figure  1.  Histograms  (black)  of  the  obtained  distributions  from  105  synthetic  data  sets  obtained  using  the  epochs  for  calibrator  2M0746+2000.  Each  simulated  data 
set  assumes  a  0  parallax  and  +100  mas  yr-1  in  both  R.A.  and  decl.  The  red  dashed  line  is  a  Gaussian  distribution  with  o  equal  to  the  standard  deviation  measured  on 
the  Monte  Carlo  generated  data  sets,  illustrating  the  perfect  agreement  between  the  two,  and  validating  our  approach  to  determine  empirical  uncertainties. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


astrometric  process  eliminating  the  following  objects  from  the 
reference  frame:  (1)  the  target  parallax  star,  (2)  other  high  proper 
motion  stars  in  the  field,  (3)  saturated  stars,  and  (4)  elongated 
or  extended  sources  (e.g.,  galaxies).  The  whole  iterative  process 
was  repeated  a  second  time.  The  final  catalog  that  was  output  by 
ATPa  contained  the  five  free  parameters  defined  above  for  each 
star  in  the  field.  It  also  contained  formal  uncertainties  derived 
from  the  covariance  matrix  of  the  least-squares  solution  for  the 
equations  above,  information  about  the  number  of  observations 
employed,  and  the  rms  of  the  residuals  per  epoch. 

In  Table  1  we  list  the  data  on  each  target  including  the  baseline 
between  the  first  and  last  observations,  the  number  of  nights  a 
target  was  observed  (epochs),  and  how  many  parallax  frames 
were  included  in  the  solution.  In  the  case  of  a  number  of  targets, 
parallactic  sampling  was  low  (e.g.,  5-7  epochs);  therefore,  to 
ensure  realistic  uncertainties,  we  applied  a  Monte  Carlo  method 
to  the  solution.  For  each  target,  we  independently  measured  the 
standard  deviation  of  the  R.A.  and  decl.  residuals.  However, 
since  many  of  the  targets  had  a  small  number  of  epochs,  we 
used  the  median  standard  deviation  from  the  reference  stars 
instead.  These  more  realistic  residuals  were  used  to  add  random 
Gaussian  noise  to  simulated  observations  with  the  parallax 
and  proper  motion  from  the  final  ATPa  catalog.  We  repeated 
this  experiment  1000  times  and  solved  for  the  five  astrometric 
parameters  each  time.  The  standard  deviation  around  the  mean 
value  of  the  parallax  over  the  Monte  Carlo  runs  is  the  final 
uncertainty  listed  in  Table  2.  As  long  as  the  standard  deviations 
used  in  the  Monte  Carlo  experiment  are  realistic,  this  approach 
automatically  accounts  for  the  correlation  between  the  parallax 
and  proper  motion  and  offsets  the  issue  of  undersampled 
astrometry.  We  flag  sources  in  Table  1  that  have  <9  epochs 
and  potentially  undersampled  parallaxes. 

We  note  that  parallax  and  proper  motion  are  linear  parameters 
of  the  model.  As  a  consequence,  one  does  not  depend  on 
the  particular  value  of  the  other  in  any  way.  Still,  it  is  true 
that  they  can  be  correlated  due  to  the  sampling  cadence.  One 
of  the  advantages  of  the  Monte  Carlo  approach  used  here  is 
that  such  a  correlation  is  automatically  accounted  for  in  the 
process  of  generating  synthetic  data  sets  at  the  same  observing 
epochs  and  deriving  the  empirical  standard  deviations  over  the 
repetitions  of  the  same  measurement.  To  illustrate  this,  we  use 
the  observing  epochs  of  2M0746+2000  (13  epochs  over  2.5  yr) 
and  simulate  astrometric  measurements  assuming  0  parallax  and 
a  proper  motion  of  +100  mas  yr-1  on  both  R.A.  and  decl.  We 
introduce  random  Gaussian  noise  with  a  standard  deviation  of 


5  mas  in  both  R.A.  and  decl.  and  solve  for  all  the  astrometric 
parameters.  We  repeat  this  experiment  10,000  times  and  produce 
the  numerical  distributions  of  the  obtained  /jLa,  /+$,  and  parallax. 
As  expected  from  a  linear  model,  the  obtained  distributions  are 
non-biased  and  the  shape  of  the  marginalized  distributions  is 
Gaussian  with  the  same  o  we  measure  on  the  scatter  of  the  105 
Monte  Carlo  obtained  parameters  (see  Figure  1). 

3.3.  Correction  from  Relative  to  Absolute  Parallax 

The  final  parallaxes  from  the  astrometric  solution  are  relative 
to  the  motion  of  the  background  stars  chosen  as  references. 
A  correction  is  required  based  upon  the  true  parallaxes  of  the 
reference  stars  to  convert  to  an  absolute  measurement. 

Typically  there  are  three  ways  to  convert  from  relative  to 
absolute  parallaxes:  (1)  using  statistical  methods  that  rely  on 
a  well-determined  model  of  the  Galaxy  and  is  most  relevant 
for  faint  distant  reference  stars,  (2)  spectroscopic  parallaxes 
that  rely  on  spectral  data  obtained  for  every  reference  star,  or 
(3)  photometric  parallaxes  that  rely  on  colors  for  all  reference 
stars.  We  determined  the  parallax  corrections  using  the  third 
method  because  the  reference  stars  are  primarily  the  brightest 
in  the  field  and  spectral  data  are  not  available. 

In  order  to  measure  photometric  parallaxes  for  the  reference 
stars  we  assume  that  all  sources  are  main-sequence  dwarfs.  Fol¬ 
lowing  a  similar  procedure  described  in  Vrba  et  al.  (2004),  we 
obtained  2MASS  photometry  for  all  reference  stars.  We  com¬ 
pared  with  the  intrinsic  colors  described  in  Koornneef  (1983) 
so  we  first  converted  J,H,KS  values  to  the  Koornneef  photo¬ 
metric  system  using  the  transformations  detailed  in  Carpenter 
(2001).  The  near-IR  J  —  H  and  H—K  colors  were  used  to  es¬ 
timate  spectral  types  and  V  —  K  colors  of  the  background  stars 
based  on  the  relations  detailed  in  Koornneef  (1983)  for  main- 
sequence  dwarfs.  Absolute  V magnitudes  were  taken  from  color 
spectral  type  relations  described  in  Kitchin  (2004)  and  then  con¬ 
verted  to  Mj  and  MK  values.  Distances  to  reference  stars  were 
determined  by  averaging  ( m  —  M)  j  and  ( m  —  M)K  values  which 
were  in  good  agreement. 

Each  reference  star  was  given  equal  weight  in  the  astrometric 
solution.  As  a  result  we  averaged  the  photometric  parallaxes  to 
calculate  the  distance  correction  and  used  a  standard  deviation 
of  the  mean  for  the  correction  uncertainty.  We  added  the  distance 
correction  to  our  relative  parallax  and  added  the  correction 
uncertainty  in  quadrature  with  our  parallax  uncertainty  to  obtain 
the  final  absolute  parallax.  The  average  correction  to  absolute 
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Table  2 

Astrometry  of  Targets 


Name 

(1) 

SpT 

(2) 

^rel 

(mas) 

(3) 

^abs 

(mas) 

(4) 

(mas  yr-1) 

(5) 

1^8 

(mas  yr-1) 

(6) 

Vtan 

(km  s_1) 
(7) 

2MASSJ00325584— 4405058 

L0  y 

36.6  ±  4.6 

38.4  ±4.8 

117.8  ±4.3 

-91.6  ±4.3 

18.4  ±2.8 

2M  ASS  J00345 1 57+0523050 

T6.5 

103.3  ±  7.2 

105.4  ±  7.5 

658.8  ±  6.2 

181.0  ±6.0 

30.8  ±  1.2 

2MASSIJ0103320+193536 

L6/3 

45.2  ±  7.5 

46.9  ±  7.6 

293.0  ±  4.6 

27.7  ±  4.7 

29.8  ±  4.2 

2MASSJ02212859— 6831400 

M8£ 

23.7  ±  3.3 

25.4  ±  3.6 

53.9  ±  4.4 

13.7  ±4.5 

10.4  ±  2.9 

2MASSJ03 185403—3421292 

L7 

70.2  ±  7.3 

72.9  ±  7.7 

391.7  ±6.4 

34.1  ±6.5 

25.6  ±  2.0 

SDSS  J032553 . 17+042540. 1 

T5.5 

52.9  ±  10.8 

55.6  ±  10.9 

-163.7  ±5.8 

-59.6  ±  5.7 

14.9  ±  3.3 

2MASSJ0334121 8—4953322 

M9 

119.1  ±3.4 

120.6  ±  3.6 

2329.5  ±  4.2 

553.6  ±  4.3 

94.2  ±  0.8 

2MASSJ04221413+1530525: 

M6  y 

22.2  ±  2.9 

24.8  ±3.1 

-17.2  ±2.7 

7.4  ±  2.6 

3.6  ±  1.5 

2MASSJ04390101 -2353083 

L6.5 

108.8  ±  3.8 

110.4  ±4.0 

-116.3  ±3.8 

-162.0  ±3.8 

8.6  ±  0.3 

2MASSJ04455387— 3048204 

L2 

75.9  ±  2.7 

78.5  ±  4.9 

164.0  ±  2.8 

-415.0  ±2.7 

27.0  ±  1.2 

2M  AS  S  J05 0 1 2406 — 00 1 045 2 

L4y 

75.2  ±  4.7 

76.4  ±  4.8 

182.4  ±4.3 

-132.7  ±4.2 

14.0  ±  0.9 

2MASSJ05160945— 0445499 

T5.5 

42.5  ±  6.4 

44.5  ±  6.5 

-202.7  ±  3.6 

-190.1  ±3.6 

29.6  ±  3.9 

2MASSJ05 184616—2756457 

LOy 

19.8  ±6.8 

21.4  ±6.9 

28.6  ±  4.2 

-16.0  ±4.0 

7.3  ±  6.6 

2MASSJ05 185995—2828372 

L7.5 

46.2  ±  6.8 

47.5  ±  6.8 

-76.5  ±  4.0 

-272.7  ±4.1 

28.3  ±  3.6 

2MASSJ05361998- 1920396 

Ll/3 

24.0  ±  9.3 

25.6  ±  9.4 

24.6  ±  5.3 

-30.6  ±  5.0 

7.3  ±7.1 

LHS1777 

M5 

77.4  ±2.1 

79.9  ±  2.5 

-175.6  ±2.2 

937.9  ±  2.3 

56.6  ±  0.8 

2MASSJ06085283— 2753583 

M8.5y 

30.9  ±  3.5 

32.0  ±  3.6 

8.9  ±  3.5 

10.7  ±  3.5 

2.1  ±0.9 

2M  AS  S  J06 1 64006 — 6407 1 94 

sdL5 

18.1  ±6.3 

19.9  ±  6.5 

1304.0  ±  3.9 

-31.9  ±3.7 

311.7  ±45.5 

2MASSJ06244595— 4521548 

L5 

83.2  ±  4.5 

83.9  ±4.5 

-41.0  ±5.4 

355.3  ±5.1 

20.2  ±  0.8 

DENIS-PJ065248. 5— 574137 

M8£ 

30.1  ±3.1 

31.3  ±3.2 

0.1  ±3.4 

29.2  ±  3.3 

4.4  ±  1.2 

2MASSJ07123786— 6155528 

Ll/8 

21.2  ±9.0 

22.9  ±9.1 

-35.7  ±  4.9 

79.1  ±4.8 

17.9  ±  13.2 

2MASSJ07290002— 3954043 

T8 

123.4  ±8.1 

126.3  ±  8.3 

-566.6  ±  5.3 

1643.4  ±  5.5 

65.3  ±  1.5 

SDSSJ074201. 41+205520.5 

T5 

65.1  ±  8.5 

66.5  ±  8.6 

-311.5  ±8.7 

-227.9  ±  8.7 

27.6  ±  2.7 

2MASSIJ0746425+200032 

L0.5 

84.8  ±  4.5 

86.2  ±  4.6 

-355.9  ±5.1 

-63.7  ±  5.2 

19.9  ±0.8 

SDSSJ083048. 80+01283 1 . 1 

T4.5 

41.2  ±5.9 

43.1  ±6.1 

188.9  ±3.0 

-351.5  ±2.9 

43.9  ±  4.6 

2MASSJ08472872- 1532372 

L2 

75.4  ±  3.5 

76.5  ±  3.5 

131.3  ±3.9 

-199.0  ±3.9 

14.8  ±  0.6 

2MASSIJ0859254- 194926 

L7 

64.2  ±  6.1 

65.4  ±6.1 

-324.5  ±  3.7 

-94.7  ±  3.5 

24.5  ±  1.8 

2MASSJ09393548— 2448279 

T8 

194.6  ±  10.4 

196.0  ±  10.4 

558.1  ±5.8 

-1030.5  ±5.6 

28.4  ±  0.6 

2MASSJ09490860- 1545485 

T2 

53.2  ±  6.4 

55.3  ±  6.6 

-98.1  ±2.9 

8.9  ±  2.9 

8.4  ±  1.5 

WT248 

M3 

30.6  ±  4.6 

31.7  ±4.7 

-1199.4  ±5.9 

-140.5  ±5.8 

180.6  ±  11.3 

2MASSJ10073369— 4555147 

T5 

69.6  ±  5.1 

71.0  ±5.2 

-723.5  ±  3.4 

148.7  ±  3.6 

49.3  ±  1.9 

2MASSIJ1010148— 040649 

L7 

58.5  ±  8.1 

59.8  ±8.1 

-308.7  ±  4.9 

-9.7  ±4.8 

24.5  ±  2.8 

2MASSJ1 0220489+0200477 

M9£ 

25.0  ±  11.4 

26.4  ±  11.5 

-156.2  ±6.6 

-429.0  ±  6.8 

82.1  ±29.8 

SDSSJ103026.78+021 306.4 

L9.5 

23.6  ±  10.6 

24.8  ±  10.6 

50.0  ±  5.3 

62.0  ±  4.8 

15.3  ±  13.1 

2MASSWJ1036530— 344138 

L6 

59.8  ±  8.9 

61.5  ±9.1 

-13.3  ±6.7 

-486.6  ±  5.9 

37.5  ±  3.7 

SDSSJ104335. 08+121 3 14.1 

L7 

66.8  ±  10.6 

68.5  ±  10.6 

26.0  ±5.1 

-234.2  ±  3.9 

16.3  ±2.5 

SDSSJ1 04409.43+042937.6 

L7 

50.0  ±  10.0 

51.3  ±  10.1 

-29.6  ±  4.3 

89.7  ±  3.4 

8.7  ±  2.7 

2MASSJ10584787- 1548172 

L3 

64.2  ±  3.9 

66.5  ±  4.4 

-253.5  ±  3.4 

39.7  ±  3.2 

18.3  ±  1.1 

TWA28 

M8.5y 

26.4  ±  4.2 

28.0  ±  4.3 

-55.7  ±  3.9 

-7.4  ±3.8 

9.5  ±  2.9 

2MASSJ1 1145133-2618235 

T7.5 

175.1  ±7.0 

176.8  ±  7.0 

-2927.2  ±  7.0 

-374.2  ±  7.2 

79.2  ±  0.8 

TWA26 

M9y 

33.6  ±  4.2 

35.1  ±4.3 

-75.2  ±  4.4 

-10.4  ±4.5 

10.3  ±2.1 

2MASSJ1 1553952-3727350 

L2 

102.8  ±  4.6 

104.4  ±  4.7 

66.8  ±4.8 

-777.9  ±  4.7 

35.5  ±  0.8 

SDSSJ1 15553.86+055957.5 

L7.5 

56.6  ±  10.2 

57.9  ±  10.2 

-420.5  ±6.1 

-54.5  ±  4.8 

34.7  ±  4.5 

SDSSJ1 15700.50+061 105.2 

T1.5 

31.4  ±  14.2 

33.8  ±  14.2 

158.9  ±7.1 

-121.3  ±6.5 

28.0  ±  14.7 

SDSSJ120747. 17+024424. 8 

L8 

42.3  ±  12.2 

44.5  ±  12.2 

-487.0  ±  12.9 

140.1  ±  11.7 

54.1  ±  10.4 

2MASSJ12095613- 1004008 

T3 

45.9  ±  11.1 

47.3  ±  11.1 

268.8  ±  7.0 

-356.7  ±  5.7 

44.8  ±  7.7 

2MASSJ12154432— 3420591 

T4.5 

38.5  ±  8.9 

39.8  ±  8.9 

-223.7  ±  5.6 

-344.8  ±  5.2 

48.9  ±  8.3 

2MASSJ135955 10— 4034582 

LI 

62.8  ±  5.4 

64.2  ±  5.5 

64.6  ±  7.7 

-494.1  ±6.7 

36.8  ±2.1 

2M  AS  S  J 1 404494 1—3159329 

T2.5 

38.5  ±6.1 

40.0  ±  6.2 

334.5  ±  4.9 

-15.5  ±4.7 

39.7  ±  5.0 

2M  AS  S  J 1 4442067 — 20 1 9222 

sdM9 

59.7  ±  4.9 

61.2  ±5.1 

-2905.1  ±3.4 

-1953.9  ±3.4 

271.3  ±5.5 

SDSSpJ14460060+002452.0 

L6 

45.2  ±  14.5 

46.7  ±  14.6 

202.9  ±  10.9 

-100.4  ±  10.8 

23.0  ±  7.7 

SDSSJ15041 1.63+102718.4 

T7 

51.2  ±7.1 

52.5  ±7.1 

373.8  ±  7.9 

-322.5  ±  7.7 

44.6  ±  4.0 

SDSSJ151 1 14.66+060742.9 

TO 

35.4  ±  6.4 

36.7  ±  6.4 

-255.6  ±7.1 

-238.0  ±7.5 

45.1  ±6.3 

SDSSJ152103. 24+013 142.7 

T2 

40.2  ±  7.2 

41.3  ±7.2 

-174.3  ±7.7 

83.7  ±  7.7 

22.2  ±  4.2 

2MASSIJ1526140+204341 

L7 

47.1  ±  8.6 

48.5  ±  8.7 

-206.6  ±  10.1 

-349.8  ±  9.0 

39.7  ±  5.3 

2MASSJ16150413+1340079 

T6 

67.2  ±  6.3 

68.6  ±  6.4 

257.8  ±  7.9 

-321.9  ±7.9 

28.5  ±  1.9 

SDSSJ163022.92+081 822.0 

T5.5 

45.4  ±  10.1 

46.9  ±  10.1 

-75.5  ±  12.5 

-100.0  ±  12.1 

12.7  ±  3.7 

2MASSJ1645221 1  —  1319516 

LI. 5 

107.6  ±  5.7 

109.9  ±6.1 

-347.7  ±  11.5 

-777.2  ±  10.4 

36.8  ±  1.0 

2MASSJ17545447+1649196 

T5.5 

86.0  ±  10.2 

87.6  ±  10.2 

113.5  ±9.1 

-141.4  ±9.2 

9.8  ±  1.3 

2MASSJ1 8283572-4849046 

T5.5 

81.0  ±7.1 

83.7  ±  7.7 

231.4  ±  10.5 

52.4  ±  10.9 

13.4  ±  1.2 

2MASSJ193601 87—5502322 

L5 

64.9  ±  5.4 

66.3  ±  5.4 

234.7  ±  9.0 

-289.3  ±  9.3 

26.6  ±  1.7 

SDSSJ204749.61— 071818.3 

TO 

48.5  ±  7.8 

49.9  ±  7.9 

48.7  ±  11.4 

-193.8  ±  11.2 

19.0  ±3.2 

SDSS  J205235 .3 1  - 1 60929.8 

T1 

25.2  ±  9.1 

26.8  ±9.1 

410.4  ±  10.9 

170.4  ±  9.7 

78.7  ±  20.9 

2MASSJ20575409— 0252302 

LI. 5 

68.4  ±  3.6 

70.1  ±3.7 

1.6  ±3.8 

-86.3  ±  3.9 

5.8  ±0.5 
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Table  2 

(Continued) 


Faherty  et  al. 


Name 

SpT 

TTrel 

^abs 

na 

ns 

Vtan 

(mas) 

(mas) 

(mas  yr-1) 

(mas  yr-1) 

(km  s-1) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

2MASSJ21321 145+1341584 

L6 

28.3 

± 

8.0 

30.0 

± 

8.2 

16.4 

± 

11.6 

-141.0 

± 

11.2 

22.4 

db 

8.0 

2MASSJ21513839— 4853542 

T4 

48.7 

± 

6.6 

50.4 

± 

6.7 

408.8 

± 

9.9 

-196.2 

± 

9.0 

42.7 

± 

4.0 

2MASSJ22282889— 43 10262 

T6 

92.7 

± 

7.0 

94.0 

± 

7.0 

102.3 

± 

5.8 

-324.4 

± 

5.1 

17.2 

± 

1.0 

2MASS  J23224684-  3 133231 

L0£ 

57.6 

± 

5.5 

58.6 

± 

5.6 

-194.8 

± 

7.4 

-527.3 

± 

7.5 

45.5 

± 

2.7 

2MASSIJ2356547— 155310 

T5.5 

71.7 

± 

4.4 

74.4 

± 

5.8 

-422.7 

± 

4.0 

-615.9 

± 

3.6 

47.7 

± 

2.0 

2MASSJ23594034— 7335055 

T5.5 

85.4 

± 

5.7 

86.7 

± 

5.7 

280.9 

± 

5.4 

39.2 

± 

5.5 

15.5 

± 

0.9 

parallax  for  the  full  list  of  targets  was  1 .5  zb  0.5  mas  ranging  from 
0.8  to  2.9  mas.  The  final  parallaxes  with  absolute  corrections 
are  shown  in  Table  2. 

3.4.  Comparison  of  Calibrators 

There  are  nine  calibrator  stars  in  our  full  astrometric  sample 
(three  imaged  with  ISPI  and  six  imaged  with  ANDICAM)  that 
we  obtained  as  a  check  on  the  reliability  of  our  methods.  Table  3 
lists  the  astrometry  for  the  calibrators  measured  in  this  work  and 
compares  those  values  with  results  reported  in  the  literature  (see 
Figure  2).  Of  the  nine  calibrators,  seven  match  within  la  and 
all  but  one  (TWA  28)  match  within  2cr  of  published  values.  No 
systematic  trends  were  detected  in  the  parallax  measurement 
differences.  The  mean  difference  and  scatter  between  literature 
and  BDKP  values  is  —4.0  ±5.0  mas.  This  indicates  that  our 
parallax  pipeline  produces  reliable  results. 

As  a  check  on  non-calibrator  stars  we  examined  how  well 
the  proper  motion  values  match  published  values  for  all  sources 
in  our  sample.  The  two  lower  panels  in  Figure  2  show  plots 
of  na  and  /x«5  from  the  literature  versus  our  calculated  values 
with  uncertainties.  We  find  that  42%  of  the  sample  match 
both  components  within  lcr,  76%  match  within  2cr,  87% 
match  within  3a,  and  all  but  3%  match  within  4a.  This  is 
different  from  a  Gaussian  error  distribution,  possibly  indicating 
underestimated  uncertainties  in  either  published  values  or  our 
own  astrometric  solution.  The  mean  difference  between  /ia  and 
fis  values  was  4  ±  32  mas  and  —4  ±  30  mas,  respectively.  No 
systematic  trends  were  detected  in  the  proper  motion  component 
differences.  We  list  objects  whose  p,a  or  fis  components  were 
discrepant  by  more  than  4a  in  Table  4.  2M1404-3159  is  among 
the  discrepant  proper  motion  objects  yet  we  believe  there  was  a 
sign  reporting  error  in  Looper  et  al.  (2007)  as  we  calculate  the 
same  magnitude  of  motion  in  the  opposite  direction. 

4.  ABSOLUTE  MAGNITUDE  RELATIONS 

As  of  2011  September  there  were  106  L  and  T  dwarfs  with 
published  parallax  measurements.  We  have  added  59  to  this 
list,  doubling  the  number  of  measurements  in  some  spectral 
bins.  The  precision  on  parallaxes  in  this  work  as  well  as  within 
the  literature  varies  greatly.  To  ensure  that  the  analysis  that 
follows  was  not  biased  by  uncertain  parallaxes  or  photometry, 
we  required  all  sources  to  have  MJHK  uncertainty  <0.5  mag.  We 
list  photometry,  parallax  measurements,  and  references  for  all 
known  and  new  UCDs  in  Table  5.  With  a  substantial  increase  in 
the  number  of  objects,  we  can  re-evaluate  the  color-magnitude 
and  spectral  type  luminosity  trends  originally  defined  by  Dahn 
et  al.  (2002),  Tinney  et  al.  (2003),  and  Vrba  et  al.  (2004), 
particularly  across  the  poorly  sampled  L/T  transition  region.  We 
list  all  L  and  T  dwarfs  with  parallaxes  and  their  corresponding 


-400  -200  0  200  400  600  800 

Literature  pa 


Figure  2.  Top  plot:  the  parallax  measurement  comparison  of  the  nine  calibrator 
dwarfs.  All  but  two  match  previous  measurements  within  lcr.  Middle  and  lower 
plots:  comparison  of  literature  proper  motion  components  to  those  measured  in 
this  work.  Three  objects  are  discrepant  by  more  than  4cr  and  listed  in  Table  4. 
In  each  plot  blue  filled  circles  represent  dwarfs  that  were  measured  with  ISPI 
and  red  five  point  stars  were  measured  with  ANDICAM. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Table  3 

Astrometric  Calibrators 


Name 

^abs 

na 

ns 

^abs 

na 

ns 

Inst.a 

Ref 

(mas) 

(mas  yr-1) 

(mas  yr-1) 

(mas) 

(mas  yr-1) 

(mas  yr-1) 

(BDKP) 

(BDKP) 

(BDKP) 

(Lit) 

(Lit) 

(Lit) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(J) 

(8) 

(9) 

LHS  1777 

79.9  ±  2.5 

-175.6  +  2.2 

937.9  ±  2.3 

78.20  ±  2.70 

-153.00  +  5.00 

957.00  ±  5.00 

A 

1 

79.72  ±  1.89 

2 

80.11  +  2.34 

3 

WT  248 

31.7  +  4.7 

-1199.4  +  5.9 

—  140.5  ±  5.8 

38.44  ±  2.83 

-1191.00  ±  13.00 

-115.00  ±  13.00 

A 

1 

2MASS  J0746+2000 

86.2  ±  4.6 

-355.9  +  5.1 

-63.7  ±  5.2 

81.90  +  0.30 

-374.04  +  0.31 

-57.91  ±  0.65 

A 

4 

2MASS  J0939— 2448 

196.0  ±  10.4 

558.1  ±5.8 

—  1030.5  ±  5.6 

187.30  +  4.60 

573.40  ±  2.30 

-1044.70  ±2.50 

I 

6 

TWA28 

28.0  +  4.3 

-55.7  ±  3.9 

-7.4  ±3.8 

18.1+0.5 

-67.2  ±  0.6 

-14.0  +  0.6 

A 

8 

2MASS  J1058— 1548 

66.5  ±  4.4 

-253.5  ±  3.4 

39.7  ±  3.2 

57.70+  1.00 

-252.93  ±  0.50 

41.42  ±0.45 

A 

4 

2MASS  J1444— 2019 

61.2  +  5.1 

-2905.1  +  3.4 

-1953.9  +  3.4 

61.67  +  2.12 

-2906.15  +  2.41 

-1963.12  +  2.71 

A 

7 

SDSS  J 1446+0024 

46.7  ±  14.6 

202.9  ±  10.9 

-100.4+10.8 

45.46  ±  3.25 

179.60  ±  6.68 

-65.58  ±  4.07 

I 

5 

2MASS  J2356— 1553 

74.4  ±  5.8 

-422.7  ±  4.0 

-615.9  +  3.6 

68.97  ±  3.42 

-443.44  ±  2.07 

-600.15  +  2.48 

I 

5 

Notes.  a  A:  ANDICAM  and  I:  ISPI. 

References.  (1)  van  Altena  et  al.  1995;  (2)  Jao  et  al.  2005;  (3)  Costa  et  al.  2005;  (4)  Dahn  et  al.  2002;  (5)  Vrba  et  al.  2004;  (6)  Burgasser  et  al.  2008c;  (7)  Schilbach 
et  al.  2009;  (8)  Teixeira  et  al.  2008. 


Table  4 

Discrepant  Proper  Motion  Values 


Name 

na 

(mas  yr-1) 

(BDKP) 

ns 

(mas  yr-1) 

(BDKP) 

na 

(mas  yr-1) 

(Lit) 

ns 

(mas  yr-1) 

(Lit) 

Reference 

2M  AS  S  J06 1 64006 — 6407 1 94 

1304.0  ±  3.9 

-31.9  +  3.7 

1405  ±  8 

-51  ±  18 

1 

2MASSJ14044941 — 3 159329a 

334.5  ±  4.9 

-15.5  +  4.7 

-348  ±  30 

32  ±3 

2 

2MASSJ21321 145+1341584 

16.4+  11.6 

-141.0+  11.2 

-55  +  9 

-395  ±  9 

3 

Notes.  Details  on  the  discrepant  proper  motion  objects.  We  note  only  objects  whose  proper  motion  values  were  discrepant  by  more  than  4a. 
a  We  find  the  same  magnitude  of  motion  as  Looper  et  al.  (2007)  yet  an  opposite  direction.  We  conclude  there  must  be  a  sign  error  in  that  work. 
References.  (1)  Cushing  et  al.  2009;  (2)  Looper  et  al.  2007;  (3)  Siegler  et  al.  2007. 


magnitudes  in  optical  through  mid-IR  wavelengths  in  Table  5. 
The  apparent  MKO  JHK  magnitudes  listed  in  Table  5  were 
used  to  calculate  the  absolute  MKO  MJHK  magnitudes  used 
throughout  the  analysis. 

Figure  3  shows  the  MJHk  magnitude  versus  spectral  type  for 
all  late-type  M,  L,  and  T  dwarfs  with  parallax  measurements. 
Seemingly  normal  field  objects  (hereafter  “normal”  is  defined  as 
excluding  tight  binaries  unresolved  in  2MASS,  young  sources 
with  low  surface  gravity  features,  and  subdwarfs)  provide  a 
guideline  for  how  near-IR  intrinsic  brightness  changes  with 
spectral  type.  Absolute  magnitudes  of  brown  dwarf  subgroups 
(unresolved  binaries  in  red,  subdwarfs  in  blue,  and  low  surface 
gravity  or  companions  to  young  field  stars  in  green  and  purple) 
allow  us  to  investigate  how  secondary  parameters  such  as 
binarity,  metallicity,  and  age  influence  brown  dwarf  observables. 
In  this  section,  we  investigate  what  can  be  extrapolated  for  the 
population  using  the  sequence  of  normal  field  dwarfs  including 
resolved  L/T  transition  binaries.  In  Sections  6,  7,  and  9  we 
discuss  in  detail  the  differences  among  the  subdwarfs,  low 
surface  gravity,  and  binary  dwarfs,  respectively. 

Brown  dwarfs  have  highly  structured  spectral  energy  dis¬ 
tributions,  and  magnitudes  in  JHK  are  extremely  sensitive  to 
the  exact  filter  bandpass  used.  Therefore,  we  converted  all 
magnitudes  onto  the  Mauna  Kea  Observatory  filter  set  (MKO; 
Tokunaga  et  al.  2002),  whose  narrow  bandpasses  are  less  af¬ 
fected  by  atmospheric  absorption  than  the  CIT  and  2MASS  fil¬ 


ter  sets  (particularly  at  J).  If  required,  the  transformations  from 
Stephens  &  Leggett  (2004)  were  used  to  convert  from  2MASS 
to  MKO  magnitudes. 

Most  of  the  L  dwarfs  in  our  sample  were  classified  spectrally 
from  red  optical  data  following  the  scheme  of  Kirkpatrick 
et  al.  (1999),  while  the  T  dwarfs  were  classified  in  the  near- 
IR  (Burgasser  et  al.  2006a).  An  optical  spectral  type  was  used 
for  any  object  classified  as  an  L  dwarf  and  a  near-IR  spectral 
type  was  used  for  any  object  classified  as  a  T  dwarf.  For  any  L 
dwarf  lacking  optical  data  we  used  its  near-IR  spectral  type. 

4.1.  Brown  Dwarf  H-R  Diagram 

Figure  4  shows  the  MJHK  sequence  for  normal  L  and  T 
dwarfs  with  uncertainties  <0.5  mag.  We  have  distinguished 
published  parallaxes  (open  circles)  from  those  reported  in  this 
work  (filled  gray  circles)  to  demonstrate  the  impact  of  these 
new  measurements.  We  have  augmented  the  field  dwarf  sample 
in  Figure  5  with  resolved  photometry  for  9  L-dwarf/T-dwarf 
transition  binaries  for  a  more  detailed  look  at  brightness 
trends  across  this  transition  (Burgasser  et  al.  2010a;  Looper 
et  al.  2008a;  McCaughrean  et  al.  2004;  Liu  et  al.  2010; 
Stumpf  et  al.  2011).  Four  known  binary  sources,  2MASS 
J0518— 2828,  2MASS  J1209-1004,  2MASS  J1404-3159,  and 
SDSS  J2052— 1609  had  parallax  measurements  reported  in  this 
work  as  did  one  suspected  binary  SDSS  J15 1 1+0607  (Burgasser 
et  al.  2010a;  C.  R.  Gelino  et  al.  in  preparation).  We  list  the 
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Table  5 

Photometry  of  Full  Sample* 


Faherty  et  al. 


Name 

(1) 

SpT 

(OpT) 

(2) 

SpT 

(near-IR) 

(3) 

7T 

(mas) 

(4) 

J 

(MKO) 

(5) 

H 

(MKO) 

(6) 

K 

(MKO) 

(7) 

Noteb 

(8) 

Ref 

(9) 

2MASSJ00043484— 4044058 

L5.0 

L4.5 

76.86  ±  3.97 

13.06  ±  0.03 

12.14  ±0.03 

11.36  ±0.03 

B 

23,  65,  9,  70 

2M  AS  S  J0025 0365 +47 59191 

L4.0 

- 

23.69  ±  1.08 

14.72  ±  0.04 

13.71  ±0.03 

12.85  ±  0.06 

B 

13,94 

2MASSJ00303013- 1450333 

L7.0 

- 

37.42  ±  4.50 

16.39  ±  0.03 

15.37  ±  0.03 

14.49  ±  0.03 

N 

15,74 

2MASSJ00325584— 4405058 

LO.Oy 

- 

38.42  ±  4.80 

14.68  ±  0.04 

13.89  ±0.03 

13.23  ±  0.04 

LG 

1,23,22 

2M  ASS  J00325937+ 1410371 

- 

L8.0 

30.14  ±5.16 

16.58  ±  0.05 

15.66  ±  0.05 

14.99  ±  0.05 

N 

34,  74 

ULASJ00340277— 0052067 

- 

T8.5 

79.60  ±  3.80 

18.15  ±0.04 

18.49  ±0.04 

18.48  ±0.04 

N 

59 

2M  ASS  J00345 1 57+0523050 

- 

T6.5 

105.36  ±  7.48 

15.11  ±0.03 

15.55  ±  0.03 

15.96  ±  0.03 

N 

1,  24,  19 

2MASSJ00361617+1821104 

L3.5 

L4.0 

114.20  ±0.80 

12.30  ±  0.03 

11.64  ±0.03 

11.04  ±0.03 

N 

10,  15,  9,  69 

HD3651B 

- 

T7.5 

90.03  ±  0.72 

16.31  ±0.03 

16.72  ±  0.03 

16.86  ±  0.03 

N 

54,  46,71 

CFBDSJ00591090— 01 14013 

- 

T9.0 

108.20  ±  5.00 

18.06  ±0.05 

18.27  ±0.05 

18.71  ±0.05 

N 

60 

Notes.  Table  contains  all  sources  with  parallaxes  reported  in  this  work  as  well  as  all  sources  used  in  the  astrometric  analysis  (with  absolute  magnitude  uncertainty 
<0.5  mag). 

a  See  the  online  version  for  the  full  table  containing  rizYJHKLM  and  IRAC  [3. 6, 4.5, 5. 8, 8.0]  bands  where  available. 

b  LG  is  a  low  surface  gravity  dwarf,  Cal  is  a  calibrator  ultracool  dwarf,  SD  is  an  ultracool  subdwarf,  and  B  is  a  tight  binary  unresolved  in  2MASS. 
c  Objects  PS043,  ULAS1416,  SOSS  J1416,  and  CFBOS  1458AB  were  added  during  the  proof  process  and  not  used  in  the  analysis  throughout  the  paper. 

References.  (1)  This  paper;  (2)  Chiu  et  al.  2006;  (3)  Cruz  et  al.  2003;  (4)  Gizis  et  al.  2000;  (5)  Golimowski  et  al.  2004;  (6)  Hawley  et  al.  2002;  (7)  Kendall  et  al.  2004; 
(8)  Kirkpatrick  et  al.  1997;  (9)  Knapp  et  al.  2004;  (10)  Reid  et  al.  2000;  (1 1)  Ruiz  et  al.  1997;  (12)  Wilson  et  al.  2001;  (13)  Cruz  et  al.  2007;  (14)  Kirkpatrick  et  al.  1999; 
(15)  Kirkpatrick  et  al.  2000;  (16)  Kirkpatrick  et  al.  2008;  (17)  Reid  et  al.  2008;  (18)  Forveille  et  al.  2004;  (19)  Burgasser  et  al.  2006a;  (20)  Burgasser  et  al.  2008a; 
(21)  Stumpf  et  al.  2008;  (22)  Cruz  et  al.  2009;  (23)  EROS  Collaboration  et  al.  1999;  (24)  Burgasser  et  al.  2003d;  (25)  Looper  et  al.  2007;  (26)  Burgasser  et  al.  2004; 
(27)  Cruz  et  al.  2004;  (28)  Burgasser  et  al.  2000b;  (29)  Burgasser  et  al.  2003b;  (30)  Cushing  et  al.  2009;  (31)  Tinney  et  al.  2005;  (32)  Burgasser  et  al.  2002a;  (33) 
Delfosse  et  al.  1997;  (34)  Geballe  et  al.  2002;  (35)  Looper  et  al.  2008a;  (36)  Burgasser  et  al.  2002a;  (37)  Lepine  et  al.  2003a;  (38)  Ellis  et  al.  2005;  (39)  Martin  et  al.  1999; 
(40)  Thorstensen  &  Kirkpatrick  2003;  (41)  Burgasser  et  al.  2004;  (42)  Fan  et  al.  2000;  (43)  Burgasser  et  al.  2003a;  (44)  Scholz  et  al.  2003;  (45)  Leggett  et  al.  2000; 
(46)  Luhman  et  al.  2007;  (47)  Lodieu  et  al.  2007;  (48)  Biller  et  al.  2006;  (49)  Strauss  et  al.  1999;  (50)  Burgasser  et  al.  1999;  (51)  Artigau  et  al.  2010;  (52)  Tsvetanov  et  al. 
2000;  (53)  Nakajima  et  al.  1995;  (54)  Mugrauer  et  al.  2006;  (55)  Pinfield  et  al.  2008;  (56)  Burgasser  et  al.  2000a;  (57)  Goldman  et  al.  2010;  (58)  Burningham  et  al.  2009; 
(59)  Warren  et  al.  2007;  (60)  Delorme  et  al.  2008;  (61)  Burningham  et  al.  2008;  (62)  Lucas  et  al.  2010;  (63)  Schneider  et  al.  2002;  (64)  Burgasser  et  al.  2007; 
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component  magnitudes  of  the  binaries  (including  all  known  bi¬ 
nary  L/T  transition  objects  with  parallaxes)  in  Table  6. 

We  have  used  the  full  parallax  sample  of  normal  objects 
with  Mjhk  uncertainties  <0.5  mag  to  re-examine  commonly 
used  near-infrared  absolute  magnitude/spectral  type  relations. 
To  determine  the  best  fit  between  these  parameters,  we  applied 
an  F- test  to  a  polynomial  with  increasing  coefficients.  The 
F- test  allows  us  to  distinguish  the  false-alarm  probability  of 
a  decreasing  x2  due  to  additional  degrees  of  freedom.  In  this 
manner,  we  converged  on  a  fourth-order  fit  to  spectral  type 
versus  MKO  MJHK  for  normal  field  dwarfs.  The  coefficients  are 
listed  in  Table  7  and  the  solution  is  overplotted  in  Figure  5. 
For  comparison,  Marocco  et  al.  (2010)  recently  reported  11 
mid-  to  late-type  T  dwarf  parallaxes  and  also  converged  upon 
a  fourth-order  polynomial  fit  to  MKO  MJHK  for  L0-T9  dwarfs. 
In  that  work,  near-IR  spectral  types  were  used  for  L0-T9  likely 
leading  to  a  steeper  fit  for  L  dwarfs  than  our  fit.  In  the  same 
manner  as  Liu  et  al.  (2006),  Marocco  et  al.  (2010)  report  two 
absolute  magnitude/spectral  type  polynomials  (overplotted  as 
red  dashed  lines  in  Figure  5):  the  first  excluding  known  binaries 
and  the  second  excluding  known  binaries  as  well  as  five  binary 
candidates  across  the  L/T  transition  (as  selected  by  Liu  et  al). 
Liu  et  al.  (2006)  found  a  0.6  mag  difference  between  the  peak  of 
the  L/T  transition  when  excluding  binary  candidates  but  with 
the  addition  of  the  Marocco  et  al.  (2010)  parallax  sample  this 


was  reduced  to  0.2  mag.  Adding  our  new  parallaxes,  we  find  no 
difference  in  a  polynomial  fit  that  excludes  or  does  not  exclude 
the  five  binary  candidates,  therefore  we  report  one  fit  for  the 
entire  sample. 

As  the  data  indicate  in  Figure  5,  the  brightening  across  the 
L/T  transition  is  more  pronounced  than  demonstrated  by  a 
polynomial  fit  from  L0  to  T9.  By  eye,  the  L  dwarfs,  L/T 
transition  objects  and  T  dwarfs  appear  to  follow  distinct  and 
independent  linear  trends.  Therefore,  as  an  alternative  to  a 
full  range  L0-T9  polynomial,  we  have  split  the  normal,  single 
source  objects  in  Figure  4  into  three  ranges  of  L0-L9,  T0-T4, 
and  T5-T9  and  fit  a  linear  polynomial  to  each.  We  obtained 
uncertainties  on  the  coefficients  by  randomly  shifting  the  known 
objects  in  each  range  within  the  given  uncertainty  10,000  times 
and  fitting  a  Gaussian  to  the  range  of  parameters.  The  best- 
fit  piecewise  lines  are  overplotted  in  Figure  4  and  reported 
in  Table  8.  For  the  L/T  transition  between  TO  and  T4,  we 
find  a  brightening  at  J  between  [1.2-1. 4]  mag,  a  brightening 
at  H  between  [0. 3-0.5]  mag,  and  a  plateau  or  dimming  at  K 
between  [—0.2  to  —0.3]  mag.  The  increased  sampling  across 
the  L/T  transition  demonstrates  that  the  brightening  is  a  real  and 
prominent  feature,  however  further  investigation  of  individual 
objects  is  necessary  to  disentangle  the  role  that  dust  and 
unresolved  binarity  play  in  creating  the  effect.  We  discuss 
the  L/T  transition  objects  in  relation  to  model  predictions  in 


10 


The  Astrophysical  Journal,  752:56  (22pp),  2012  June  10 


Faherty  et  al. 


Figure  3.  Spectral  type  vs.  absolute  magnitude  in  the  MKO  JHK  filters  for 
normal  (black  filled  circles),  combined-light  binary  (red  filled  circles),  young 
(companions  to  known  < ^  1  Gyr  stars;  purple  filled  circles),  low  surface  gravity 
(green  filled  circles),  and  subdwarf  (blue  filled  circles)  M  through  T  dwarfs. 
All  UCDs  with  a  parallax  measurement  in  this  work  or  within  the  literature  are 
shown  regardless  of  uncertainty. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Section  5.3.  We  have  used  the  full  parallax  sample  of  the 
normal  objects  with  MJHK  uncertainties  L0.5  mag  to  re-examine 
commonly  used  near-infrared  absolute  magnitude/spectral  type 
relations.  Details  on  the  full  sample  of  objects  used  in  our 
analysis  as  well  as  mean  absolute  magnitude  values  in  MJHK 
are  shown  in  Table  9 

4.2.  Color-Magnitude  Trends  for  L  and  T  Dwarfs 

We  have  collected  photometric  information  (from  optical 
to  mid-IR)  for  all  known  L  and  T  dwarfs  with  parallaxes 


Spectral  Type 

Figure  4.  Spectral  type  vs.  absolute  magnitude  in  the  MKO  JHK  filters  for 
normal  L  through  T  dwarfs.  Unfilled  circles  are  ultracool  dwarfs  with  parallax 
measurements  gathered  from  the  literature  and  filled  circles  are  those  reported 
in  this  work.  Low  surface  gravity,  subdwarf,  binary,  and  young  companion 
brown  dwarfs  are  not  shown.  Only  objects  with  absolute  magnitude  uncertainties 
<0.5  mag  are  displayed. 


(106  from  the  literature  and  59  added  in  this  work),  totaling 
165  (see  Table  5).  We  examined  various  combinations  of 
optical,  near-IR,  and  mid-IR  colors  to  find  correlations  on 
color-magnitude  diagrams  that  provided  the  strongest  insight 
into  differentiating  brown  dwarf  spectral  types  and/or  physical 
properties  of  the  population.  Figures  6  and  7  show  representative 
color-magnitude  diagrams  using  only  objects  whose  absolute 
magnitude  uncertainties  are  <0.5  mag. 

Figure  6  shows  the  most  striking  linear  relationship  between 
color  and  absolute  magnitude.  Patten  et  al.  (2006)  first  showed  a 
relatively  smooth  progression  of  M  through  T  dwarfs  in  the  M Ks 
versus  Ks-[4.5]  diagram.  We  verify  that  this  relation  separates 
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Spectral  Type 

Figure  5.  Spectral  type  vs.  absolute  magnitude  in  the  MKO  JHK  filters  for 
normal  L  and  T  dwarfs  and  L/T  transition  binaries  with  resolved  components. 
Light  gray  filled  circles  represent  normal  field  dwarfs.  Black  filled  circles  are  the 
resolved  absolute  magnitudes  of  the  binary  components.  Overplotted  on  each 
panel  is  the  best-fit  polynomial  for  L0-T9  dwarfs  (see  Table  7)  in  solid  black, 
the  bright/faint  polynomials  in  red-dashed  from  Marocco  et  al.  (2010),  and  the 
best-fit  linear  functions  to  the  ranges  L0-L9,  T0-T4,  and  T5-T9  in  solid  blue. 
Only  objects  with  absolute  magnitude  uncertainties  <0.5  mag  are  shown. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


the  early  and  mid-L  dwarfs  (0.5  <  K-[ 4.5]  <  1.4)  from  the 
early  and  mid-T  dwarfs  (1.4  <  K-[4.5]  <  4.0).  However,  it 
does  not  show  a  linear  progression  with  spectral  subtypes.  The 
degeneracy  is  most  clearly  depicted  for  the  T6.0  and  T6.5  dwarfs 
(designated  by  gray  triangle  symbols)  which  have  a  nearly 
2.0  mag  spread  in  K-[4.5]  color. 

Interestingly  the  low  surface  gravity  dwarfs,  depicted  as 
filled  downward-facing  triangles,  appear  slightly  overlumi- 
nous  for  their  color  in  comparison  to  normal  field  dwarfs. 


Figure  6.  Tightest  linear  relationship  between  color  and  luminosity  identified 
for  brown  dwarfs.  L  dwarfs  smoothly  transition  into  T  dwarfs  although  the 
figure  is  nonlinear  for  spectral  subtypes  (as  depicted  by  the  T6  dwarfs  in  gray 
upward-facing  triangles).  Unresolved  binaries  are  shown  as  open  circles.  Low 
surface  gravity  dwarfs  are  shown  as  downward-facing  triangles  and  appear 
slightly  overluminous  for  their  color. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

This  is  in  contradiction  to  findings  discussed  in  Section  6 
where  objects  are  underluminous  in  MJHK  for  their  spectral 
type.  The  discrepancy  between  the  behavior  of  low  grav¬ 
ity  objects  on  the  mid-IR  color-magnitude  diagram  and  the 
spectral-type/ absolute  magnitude  diagrams  suggests  that  spec¬ 
tral  type  is  a  poor  gauge  for  effective  temperature  of  young 
objects. 

The  most  prevalent  change  on  these  color-magnitude  dia¬ 
grams  is  a  color  reversal.  Figure  7  shows  the  three  most  studied 
versions  of  this  effect.  Condensate  cloud  opacity  is  likely  the 
dominate  contributor  to  the  increasingly  red  colors  of  L  dwarfs 
and  the  onset  of  CH4  is  responsible  for  the  change  to  blue  as 
temperatures  drop  into  the  T  dwarf  regime.  As  objects  transition 
from  L  into  T  dwarfs  the  brightening  (discussed  in  Section  4.1) 
is  most  clearly  seen  in  the  J  band  although  it  is  seen  to  a  lesser 
effect  in  H  and  plateaus  in  K.  Subdwarfs  occupy  their  own  space 
on  this  diagram  with  extremely  blue  near-IR  colors.  Low  surface 
gravity  dwarfs  fall  within  the  color  space  of  normal  field  dwarfs; 
however,  they  are  redward  and  underluminous  compared  to  nor¬ 
mal  dwarfs  with  the  same  spectral  type.  In  Section  5  below,  we 
discuss  the  MK  versus  J  —  K  color-magnitude  diagram  in  con¬ 
text  against  models  to  disentangle  subtle  effects  that  drive  much 
of  the  scatter  seen  in  Figure  7. 

5.  COMPARISON  TO  EVOLUTIONARY  MODELS 

To  put  the  observed  trends  on  color-magnitude  diagrams  in 
context,  we  compared  the  data  to  two  sets  of  evolutionary  mod¬ 
els.  Saumon  &  Marley  (2008)  present  a  set  of  models  that  in¬ 
clude  a  cloud  sedimentation  parameter  and  three  gravity  choices 
(log(g)  =  [4. 5, 5. 0,5. 5])  which  can  be  varied  to  explain  (with  dif¬ 
ferent  levels  of  accuracy)  the  near-IR  color-magnitude  diagram 
for  L  and  T  dwarfs.  Burrows  et  al.  (2006)  present  a  model 
that  includes  refractory  clouds  as  well  as  a  completely  cloud¬ 
less  model  with  varying  gravity  and  metallicity  parameters.  In 
Figures  8  and  9  we  examine  the  MK  versus  J  —  K  diagram  for 
L  and  T  dwarfs  using  the  full  sample  with  trigonometric  paral¬ 
laxes  and  the  respective  evolutionary  models.  In  Figure  8,  the 
sedimentation  parameter  (/sed)  from  Saumon  &  Marley  (2008) 
is  shown  with  increasing  value  and  varying  gravity  to  repre¬ 
sent  decreasing  cloud  thickness.  In  Figure  9,  metallicity  and 
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Table  6 

Magnitudes  for  L/T  Transition  Binaries 


Name 

(1) 

SpT 

Combined 

(2) 

SpT 

Primary 

(3) 

SpT 

Secondary 

(4) 

Mj 

Primary 

(5) 

Mj 

Secondary 

(6) 

Mh 

Primary 

(7) 

Mh 

Secondary 

(8) 

Mk 

Primary 

(9) 

Mk 

Secondary 

(10) 

Ref 

(ID 

SD0423— 0414 

L7.5 

L6.5 

T2.0 

13.95  +  0.10 

14.39  +  0.12 

13.05  +  0.10 

13.78  +  0.12 

12.38  ±  0.09 

13.51  +  0.06 

1 

2M05 18-2828 

L7.5 

L6.0 

T4.0 

14.34  ±  0.41 

14.47  ±  0.42 

12.96  ±  0.39 

14.29  +  0.51 

12.14  +  0.37 

14.32  ±  0.39 

1 

SD1021— 0304 

T3.0 

T1.0 

T5.0 

14.24  ±  0.48 

14.40  +  0.51 

13.42  ±  0.40 

14.54  ±  0.57 

13.19  +  0.36 

14.65  ±  0.35 

1 

2M1404— 3159 

T2.5 

T1.0 

T5.0 

14.47  ±  0.46 

13.94  +  0.43 

13.57  ±  0.43 

14.05  ±  0.45 

12.87  ±  0.45 

14.07  ±  0.38 

2 

2M17 11+2232 

L6.5 

L5.0 

T5.5 

15.08  ±  0.52 

16.00  ±  0.59 

13.54  ±  0.44 

15.64  ±  0.75 

12.39  ±  0.42 

15.44  ±  0.65 

1 

SD 15 11+0607 

TO.O 

L5.5 

T5.0 

14.17  +  0.48 

14.72  ±  0.55 

13.18  +  0.44 

14.96  ±  0.65 

12.43  ±  0.42 

15.09  ±  0.47 

1 

Epsilon  Indi 

T2.5 

T1.0 

T6.0 

14.69  ±  0.04 

15.63  ±  0.04 

13.99  ±  0.04 

15.75  ±  0.04 

13.51+0.04 

15.69  ±  0.01 

3 

2M 1209— 1004 

T3.0 

T2.0 

T7.5 

14.17  +  0.57 

15.64  ±  0.69 

13.69  ±  0.56 

16.45  ±  0.78 

13.59  +  0.56 

16.85  +  1.04 

4 

2M2052— 1609 

T1.0 

T1.0 

T3.5 

14.04  ±  0.92 

14.19  ±  0.93 

13.05  ±  0.89 

13.61  +  0.89 

12.67  ±  0.92 

13.42  +  0.79 

5 

References.  (1)  Burgasser  et  al.  2010a;  (2)  Looper  et  al.  2008a;  (3)  McCaughrean  et  al.  2004;  (4)  Liu  et  al.  2010;  (5)  Stumpf  et  al.  2011. 

Table  7 

Coefficients  of  Polynomial  Fits  for  L0-T8  Dwarfs 


MKO  Filter 
(1) 

co 

(2) 

ci 

(3) 

C2 

(4) 

C3 

(5) 

c4 

(6) 

rms 

(7) 

Mja 

1.48948e+l 

—  1.99967e+0 

2.75734e-l 

—  1.29143e-2 

2.03252e-4 

0.30 

Mh a 

1.32372e+l 

—  1.60036e+0 

2.27363e-l 

—  1.08205e-2 

1.75368e-4 

0.27 

MK a 

9.76100e+0 

— 4.85519e-l 

9.408 16e-2 

— 4.69032e-3 

8.16516e-5 

0.28 

Notes.  Relations  use  (MKO)  magnitudes.  Polynomial  fits  to  optical  L  dwarfs  and  NIR  T  dwarfs  (L  dwarfs  with  no  optical  spectral  type 
have  NIR  spectral  types)  excluding  subdwarfs,  low  gravity  dwarfs,  and  binaries.  Function  is  defined  as  Mj,h,k  =  Yfi= o  Q(SpT)'  and  is 
valid  for  spectral  types  L0-T8  where  10  =  L0,  20  =  TO,  etc. 

a  Including  8  L/T  transition  binaries  with  resolved  photometry:  SDSS  J0423— 0414,  2MASS  J0518— 2828,  SDSS  J102 1—0304,  2MASS 
J1404— 3159,  2MASS  J1711+2232,  SDSS  J1511+0607,  Epsilon  Indi,  SDSS  J1534+1615. 


Table  8 

Coefficients  of  Piecewise  Fits  for  L0-T9  Dwarfs 


MKO  Filter 

(1) 

co 

(2) 

ci 

(3) 

SpT  range 

Mj 

0.366  ±  0.002 

8.003  ±  0.025 

L0-L9 

Mh 

0.335  ±  0.002 

7.616  ±  0.025 

Mk 

0.297  ±  0.002 

7.500  ±  0.025 

Mj 

-0.304  ±  0.020 

21.039  +  0.438 

T0-T4.0 

Mh 

-0.089  ±  0.020 

15.797  ±  0.432 

Mk 

0.082  ±  0.019 

11.729  +  0.430 

Mj 

0.892  ±  0.008 

-8.035  ±  0.212 

T5.0-T9.0 

Mh 

0.977  ±  0.008 

-10.119  +  0.225 

Mk 

0.895  ±  0.006 

-7.632  +  0.163 

Note.  See  Table  7  comments. 


gravity  variations  are  examined  using  the  models  of  Burrows 
et  al.  (2006).  All  early  and  late-type  L  and  T  dwarfs  (including 
subdwarfs)  with  trigonometric  parallaxes  ( oMk  <0.5  mag)  are 
overplotted.  In  the  following  subsections,  we  discuss  how  vari¬ 
ations  of  the  models  fit  different  spectral  types.  The  low  surface 
gravity  dwarfs  are  discussed  in  detail  in  Section  6. 

5.1.  L  Dwarfs 

Varying  gravity  and  metallicity  within  the  cloud  model  of 
Burrows  et  al.  (2006)  encompasses  the  majority  of  early  L 
dwarfs  (top  right  and  bottom  two  panels  of  Figure  9);  how¬ 
ever,  late-type  L  dwarfs  are  still  poorly  represented.  Compared 
to  the  highest  gravity,  super  solar  metallicity  track,  there  are  a 
number  of  late-type  L  dwarfs  that  are  fainter  and  redward  of 
predictions.  In  Figure  8,  the  L  dwarf  sequence  is  best  modeled 
with  the  /sed  =  1,2  parameters  (top  right  and  bottom  left  plots 


of  Figure  8).  However,  there  are  a  handful  of  red  or  potentially 
“ultra-cloudy”  objects  that  are  not  fit  by  either  model.  Significant 
outliers  include  2MASS  J 1442+6603,  which  is  a  close 
(~30  AU)  companion  to  the  moderately  young  Ml. 5  dwarf 
G239-25  (Forveille  et  al.  2004),  and  2MASS  J0619-5803 
which  is  a  companion  (~260  AU)  to  the  young  K2  star  AB 
Pic  (Chauvin  et  al.  2005).  The  independent  assessment  that 
these  objects  are  young  (or  moderately  young)  yet  redward  of 
the  cloudy  model  implies  a  connection  between  a  youth  and  a 
dusty  photosphere  (further  discussed  in  Section  6). 

5.2.  T  Dwarfs 

The  mid-  to  late-type  T  dwarfs  are  best  fit  by  the  thin 
clouds  (/Sed  =  4)  track  from  Saumon  &  Marley  (2008)  and 
the  clear  model  from  Burrows  et  al.  (2006).  In  the  case  of 
the  latter  model  (bottom  right  panel  of  Figure  8)  the  predicted 
range  in  both  MK  and  J  —  K  shows  very  little  spread  whereas 
empirical  measurements  show  significant  scatter.  There  are  a 
handful  of  T  dwarfs  including  those  reported  in  this  work, 
2MASS  J1 114-2618,  and  2MASS  J1754+1649,  as  well  as 
previously  reported  ULAS  J0034-0052,  ULAS  J1335+1130, 
CFBDS  J0059— 0114,  ULAS  J0722-0540,  and  ROSS  458C 
which  are  notably  underluminous  and  red  compared  to  the 
/sed  =  4  model  predictions.  The  colors  of  these  late-type  T 
dwarfs  are  better  (or  equally)  fit  by  the  /se d  =  2  parameter 
(thicker  clouds)  which  also  encompasses  the  majority  of  mid- 
to  late-type  L  dwarfs.  Comparing  the  spread  to  the  Burrows 
et  al.  (2006)  clear  model  (top  left  panel  of  Figure  9)  shows 
similar  red,  underluminous  outliers.  Inconsistencies  with  both 
models  suggest  that  thick  condensate  clouds  continue  to  play  a 
role  in  the  photospheres  of  some  cooler  dwarfs  (see  discussion 
in  Burgasser  et  al.  2010b;  Marley  et  al.  2010).  At  least  one  of  the 
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Table  9 

Absolute  Magnitudes  With  Resolved  L/T  Binaries 


SpT 

BDKP 

Lit 

Total 

Norm 

Bin 

LG,Y 

SD 

Mj 

Mh 

Mk 

(MKO) 

(MKO) 

(MKO) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

M6 

1 

8 

9 

8 

0 

1 

0 

10.23  ±  0.43 

9.67  ±  0.42 

9.33  ±  0.41 

M7 

0 

7 

7 

4 

0 

0 

3 

10.93  ±  0.20 

10.28  ±0.13 

9.90  ±  0.09 

M8 

4 

5 

9 

4 

0 

4 

1 

10.87  ±  0.27 

10.22  ±  0.24 

9.84  ±  0.22 

M9 

2 

5 

7 

3 

0 

1 

3 

11.72  ±0.48 

11.06  ±0.55 

10.58  ±  0.55 

L0 

2 

8 

10 

6 

2 

2 

0 

11.69  ±0.24 

11.05  ±0.21 

10.51  ±0.19 

LI 

3 

5 

8 

7 

0 

1 

0 

12.26  ±  0.34 

11.52  ±0.37 

10.92  ±  0.41 

L2 

3 

2 

5 

4 

1 

0 

0 

12.82  ±  0.03 

12.09  ±  0.07 

11.44  ±0.16 

L3 

0 

6 

6 

4 

2 

0 

0 

12.84  ±0.15 

11.93  ±0.10 

11.33  ±0.09 

L4 

1 

9 

10 

6 

2 

1 

1 

13.15  ±0.31 

12.26  ±0.31 

11.60  ±0.34 

L5 

3 

7 

10 

8 

2 

0 

0 

13.51  ±0.25 

12.62  ±  0.22 

11.99  ±0.26 

L6 

4 

4 

8 

6 

1 

1 

0 

14.52  ±  0.34 

13.42  ±  0.29 

12.47  ±  0.32 

L7 

8 

7 

15 

9 

5 

0 

1 

14.44  ±  0.32 

13.55  ±0.31 

12.78  ±  0.33 

L8 

0 

9 

9 

8 

1 

0 

0 

14.66  ±  0.35 

13.70  ±0.31 

13.00  ±  0.29 

TO 

2 

0 

2 

1 

1 

0 

0 

- 

- 

- 

T1 

2 

4 

6 

5 

1 

0 

0 

14.24  ±  0.29 

13.42  ±0.41 

13.19  ±0.47 

T2 

4 

3 

7 

5 

1 

1 

0 

14.30  ±  0.37 

13.77  ±  0.23 

13.51  ±0.25 

T3 

1 

2 

3 

2 

1 

0 

0 

14.19  ±0.18 

13.74  ±0.10 

13.82  ±  0.28 

T4 

3 

1 

4 

4 

0 

0 

0 

14.22  ±  0.52 

14.34  ±  0.52 

14.55  ±  0.47 

T5 

9 

5 

14 

13 

1 

0 

0 

14.70  ±  0.58 

14.89  ±  0.53 

14.82  ±0.51 

T6 

3 

11 

14 

13 

1 

0 

0 

15.34  ±  0.27 

15.54  ±  0.30 

15.46  ±  0.52 

T7 

2 

8 

10 

10 

0 

0 

0 

16.02  ±  0.72 

16.49  ±  0.77 

16.63  ±  0.87 

T8 

1 

4 

5 

5 

0 

0 

0 

16.53  ±  0.69 

16.91  ±0.80 

17.04  ±  1.08 

T9 

0 

3 

3 

3 

0 

0 

0 

18.23  ±0.34 

18.44  ±0.33 

18.88  ±0.44 

Notes.  All  objects  with  parallaxes  and  absolute  mag  uncertainty  <0.5.  Note  that  Bin  are  binaries,  SD  are  subdwarfs,  LG  are  low  surface  gravity  dwarfs,  and  Y  are 
young  low  mass  companions. 


ultra-cloudy  T  dwarf  outliers,  Ross  45 8C,  is  known  to  be  young 
(Burgasser  et  al.  2010b;  Bumingham  et  al.  2011).  A  connection 
between  youth  and  dust  in  the  photospheres  of  T  dwarfs  is 
consistent  with  the  L  dwarf  results  discussed  in  Sections  5.1 
and  6  of  this  work  as  well  as  recent  modeling  work  of  young 
exoplanets  (Barman  et  al.  2011;  Currie  et  al.  2011). 

5.3.  L/T  Transition  Dwarfs 

The  L/T  transition  objects  are  not  fit  by  a  single  /se d 
parameter  using  the  Saumon  &  Marley  (2008)  models  nor  by 
any  single  combination  of  gravity  and/or  metallicity  on  the 
Burrows  et  al.  (2006)  models.16  The  steady  and  significant 
decrease  in  J  —  K  color  with  near  constant  MK  for  the  objects 
in  the  transition  region  has  been  attributed  to  the  clearing  of 
clouds  or  a  change  in  the  atmospheric  dynamical  state  (e.g., 
overall  cloud  thickness)  as  temperatures  cool  into  the  T  dwarfs 
(see  Burgasser  et  al.  2002b;  Golimowski  et  al.  2004;  Knapp 
et  al.  2004).  Without  clouds  to  provide  a  significant,  nearly  gray 
opacity,  flux  can  emerge  through  molecular  opacity  windows  in 
J  and  //bands  explaining  the  significant  brightening  discussed  in 
Section  4.1. 

To  investigate  whether  our  expanded  parallax  sample  supports 
a  rapid  cloud  clearing,  we  created  a  hybrid  model  using  the 
Saumon  &  Marley  tracks.  Similar  to  the  work  of  Burgasser  et  al. 
(2002b)  we  varied  the  sedimentation  parameter  between  the 
/sed  =  2  and  /se =  4  models  across  the  region  between  13.0< 
Mk  <  15.0.  We  started  with  the  /se d  =  2  color,  then  added  the 
/sed  =  4  color  in  10%  increments  across  the  transition.  The  result 
is  plotted  in  Figure  10.  The  L/T  transition  objects  lie  within  an 


16  Burrows  et  al.  (2006)  also  note  that  variations  in  cloud  particle  size  cannot 
account  for  the  transition  objects. 


absolute  magnitude  range  corresponding  to  model  temperatures 
spanning  ±150  K  from  a  mean  Tcu  that  depends  on  the  gravity 
chosen. 

For  an  intermediate  surface  gravity  of  log(g)  =  5.0,  which  is 
consistent  with  field  age  objects  (—3  Gyr),  the  mean  Te ff  of  our 
model  is  —  1200  ±  100  K  for  T0-T4  dwarfs.  Golimowski  et  al. 
(2004)  empirically  measured  the  Te ff  of  1 1  L7-T4  dwarfs  and 
reported  a  near  constant  Te ff  of  —  1450  K.  This  warmer  constant 
temperature  across  a  broader  range  of  spectral  types  was  strong 
evidence  for  unresolved  binarity  across  the  transition.  Indeed,  2 
of  the  1 1  objects  in  that  sample  have  since  been  resolved  into 
near-equal  mass  binaries.  The  splitting  of  spectral-type/ absolute 
magnitude  polynomials  into  bright  and  faint  samples  (see 
Section  4.1)  was  largely  driven  by  this  result.  However,  our 
result  is  in  agreement  with  the  outcome  of  other  toy  and 
sophisticated  models  (e.g.,  Burgasser  et  al.  2002b;  Saumon 
&  Marley  2008;  Marley  et  al.  2010)  which  converged  upon 
a  similar  cooler  Te ff .  Our  comparatively  lower  temperature  than 
found  by  Golimowski  et  al.  (2004)  across  the  transition  is 
consistent  with  our  finding  that  the  brown  dwarf  temperature 
plateau  occurs  across  a  narrower  spectral  subtype  range  (T0-T4 
rather  than  L7-T4)  than  previously  thought. 

6.  LOW  SURFACE  GRAVITY  DWARFS 

A  subset  of  our  parallax  sample  are  the  low  surface  gravity 
dwarfs  including  seven  M  and  seven  L  dwarfs.  Their  optical 
spectra  are  characterized  by  unusually  weak  FeH  absorption, 
weak  Na  i  and  K  i  doublets,  and  very  strong  vanadium  oxide 
bands  (Cruz  et  al.  2009;  Kirkpatrick  et  al.  2006).  They  have 
extreme  red  near-IR  colors  and  small  tangential  velocities 
relative  to  the  rest  of  the  brown  dwarf  population  (Faherty  et  al. 
2009).  Kirkpatrick  et  al.  (2006)  and  Cruz  et  al.  (2009)  have 
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Figure  7.  Near-IR  color  vs.  absolute  magnitude  diagrams  in  a  combination  of 
MKO  JHK  filters.  Spectral  subtypes  are  color-coded  and  low-gravity  as  well 
as  subdwarfs  are  depicted  as  downward-facing  triangles  and  five  point  stars, 
respectively. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


suggested  that  a  number  of  the  low  surface  gravity  dwarfs  are 
candidate  members  of  nearby  moving  groups  such  as  Pictoris, 
Tucana-Horlogium,  and  AB  Doradus,  implying  ages  roughly 
spanning  10-50  Myr  (e.g.,  Rice  et  al.  2010). 

Younger  UCDs  have  been  examined  on  color-magnitude 
diagrams,  but  as  only  a  handful  have  reported  parallaxes,  reports 
of  under-  or  -over  luminosity  have  been  speculative  and  in  some 
cases  contradictory.  There  is  evidence  that  young  dwarfs  at 
the  L/T  transition  such  as  HD  203030B  and  HN  Peg  B  are 
underluminous  or  have  a  lower  Teff  than  equivalent  spectral 
type  dwarfs.  This  has  been  explained  as  a  gravity  dependent 
temperature/ spectral  type  relation  at  the  transition  (Luhman 
et  al.  2007;  Metchev  &  Hillenbrand  2006).  HR  8799b  as  well  as 


the  early  L  dwarf  companion  to  AB  Pic  are  also  underluminous 
while  some  earlier  L  and  M  dwarfs  such  as  HD  130948B  and 
CD- 3 5  2722  B  are  overluminous  on  color-magnitude  diagrams 
(Bowler  et  al.  2010;  Chauvin  et  al.  2005;  Dupuy  et  al.  2008; 
Wahhajetal.  2011). 

Figure  11  shows  the  near-IR  absolute  magnitude  versus 
spectral  type  diagrams  for  normal  mid-type  M  through  late- 
type  L  dwarfs  with  the  low  surface  gravity  dwarfs  and  young 
companions  overplotted.  Table  10  lists  the  absolute  magnitude 
for  each  object  in  JHK  as  well  as  the  deviation  from  the  MJHK 
values  for  each  spectral  bin  calculated  from  the  polynomial 
in  Table  7.  Of  the  seven  M  dwarfs  in  this  sample,  three  are 
>0.5  mag  overluminous  for  their  spectral  type  as  might  be 
expected  for  a  young  object  which  has  not  contracted  to  its  final 
radius.  Indeed,  two  of  the  three  overluminous  M  dwarfs  are 
suspected  members  of  the  TW  Hydrae  association.  However, 
eight  out  of  the  ten  low  gravity  or  young  companion  L  dwarfs 
(~80%)  are  [0.2-1. 0]  mag  underluminous  for  the  average  of 
their  spectral  subtype  (compared  to  the  polynomial  described 
in  Table  7)  in  one  or  more  near-IR  bands.  In  Figure  11,  we 
have  plotted  the  /3  and  y  designations  assigned  for  each  object 
to  indicate  intermediate  and  low  gravity,  respectively  (see  the 
discussion  in  Kirkpatrick  2005;  Cruz  et  al.  2009).  Within  this 
sample,  there  does  not  appear  to  be  a  correlation  between  AMjhk 
(defined  as  the  difference  in  MJHK  between  the  source  and  the 
predicted  polynomial  value)  and  the  strength  of  low  surface 
gravity  features. 

The  trend  of  low  surface  gravity  UCDs  appearing  underlu¬ 
minous  for  their  spectral  type  is  surprising  given  that  young 
M  dwarfs,  such  as  those  in  TW  Hydrae,  are  1-2  mag  overlumi¬ 
nous  (D.  L.  Looper  et  al.  in  preparation).  While  we  find  that  the 
TW  Hydrae  M  dwarfs  in  our  sample  are  indeed  overluminous, 
the  L  dwarfs  show  a  different  trend.  According  to  the  evolution¬ 
ary  tracks  of  Burrows  et  al.  (1997),  10  Myr  objects  with  masses 
ranging  from  10  to  75  MJup  have  radii  that  are  25%-75%  larger 
than  1-3  Gyr  dwarfs  with  equivalent  temperatures.  This  trans¬ 
lates  into  an  over  luminosity  of  0.5-1 .2  mag.  For  50  Myr  objects 
radii  can  be  13%-50%  larger  and  would  be  0.3-0. 9  mag  over- 
luminous.  We  speculate  that  there  are  at  least  two  factors  that 
could  contribute  to  the  underluminosity:  first,  the  low-gravity 
spectral  classification  scheme  may  have  a  different  temperature 
relation  than  the  Kirkpatrick  et  al.  (1999)  classification  scheme 
used  for  normal  field  dwarfs.  For  example,  an  L0 y  dwarf  might 
have  a  significantly  different  temperature  (and  luminosity)  than 
a  normal  L0  dwarf.  Second,  young  objects  could  have  dustier 
photospheres  than  field- aged  dwarfs  thus  making  young  objects 
appear  both  fainter  and  redder  than  field  objects  of  similar  tem¬ 
perature.  Observationally,  both  low  surface  gravity  dwarfs  and 
dusty  L  dwarfs  show  red  near-IR  colors  and  similar  spectral 
characteristics  (Looper  et  al.  2008b;  Allers  et  al.  2010).  Evolu¬ 
tionary  models  demonstrate  that  the  lower  gravity  and  dustier 
(lower  /sed)  tracks  have  redder  near-IR  colors  than  intermediate, 
high  gravity,  or  larger  /se d  tracks  (see  Figure  12). 

In  Figure  12,  we  isolate  the  low  surface  gravity  L  dwarfs  with 
Mk  uncertainties  <0.5  mag  on  a  color-magnitude  diagram  with 
the  Saumon  &  Marley  (2008)  and  Burrows  et  al.  (2006)  evolu¬ 
tionary  tracks  overplotted.  In  general  the  sources  do  not  follow 
the  low-gravity  track.  Moreover,  each  model  traces  objects  at 
temperatures  that  are  several  hundred  degrees  lower  than  ex¬ 
pected  for  equivalent  spectral  type  objects.  For  example,  temper¬ 
atures  for  L4  dwarfs  range  from  ~  1600  to  1900  K  (Golimowski 
et  al.  2004);  however,  the  L4 y  dwarf,  2MASS  J0501— 0010,  is 
traced  by  model  temperatures  of  ~  1200-1 300  K. 
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Figure  8.  J  —  K  vs.  Mk  diagrams  for  L  and  T  dwarfs  with  the  evolutionary  models  of  Saumon  &  Marley  (2008)  overplotted.  Varying  the  cloud  thickness  parameter, 
/sed>  from  thick  (/se d  =  1)  to  thin  (/se d  =  4)  to  cloudless  (denoted  as  /se d  =  0),  fits  the  L  and  T  dwarf  sequences  with  varying  degrees  of  accuracy.  The  three  different 
tracks  in  each  plot  represent  different  gravities  (log(g)  =  [4. 5, 5.0, 5. 5]). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  10 

Low  Gravity  Dwarfs 


Name 

Spt 

Mj 

Mh 

Mk 

Am/ 

A  Mh* 

A  uKa 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(J) 

(8) 

2MASS  J0422+1530 

M6y 

9.6  ±  0.3 

8.8  ±0.3 

8.2  ±  0.3 

0.7 

1.0 

1.1 

2MASS  J0221— 6831 

M8£ 

10.9  ±  0.3 

10.3  ±  0.3 

9.8  ±  0.3 

-0.2 

-0.2 

0.0 

2MASS  J0652— 5741 

M8/3 

11.0  ±0.3 

10.5  ±  0.2 

9.9  ±  0.2 

-0.3 

-0.3 

-0.1 

2MASS  J0608— 2753b 

M8.5y 

11.1  ±0.3 

10.5  ±  0.3 

9.9  ±  0.3 

-0.1 

-0.1 

0.1 

TWA28 

M8.5y 

10.2  ±  0.4 

9.6  ±  0.4 

9.1  ±0.4 

0.8 

0.7 

0.9 

2MASS  J 1022+0200 

M9/3 

11.1  ±  1.0 

10.5  ±  1.0 

10.0  ±  1.0 

0.0 

0.0 

0.2 

TWA  26 

M9  y 

10.3  ±  0.3 

9.8  ±  0.3 

9.2  ±  0.3 

0.8 

0.8 

0.9 

2MASS  J05 18-2756 

L0  y 

11.8  ±0.8 

11.0  ±0.8 

10.2  ±  0.8 

-0.2 

-0.1 

0.2 

2MASS  J0032— 4405 

L0  y 

12.6  ±  0.3 

11.8  ±0.3 

11.2  ±0.3 

-1.0 

-0.9 

-0.7 

2MASS  J2322— 3133 

L0£ 

12.3  ±  0.2 

11.7  ±0.2 

11.1  ±0.2 

-0.8 

-0.8 

-0.7 

2MASS  J07 12-6155 

IAP 

12.0  ±  1.0 

11.2  ±  1.0 

10.4  ±  1.0 

0.1 

0.1 

0.3 

2MASS  J0536— 1920 

IAP 

12.7  ±  0.9 

11.8  ±0.9 

10.9  ±  0.9 

-0.6 

-0.5 

-0.1 

AB  Picb 

Llc 

12.8  ±0.2 

11.4  ±0.2 

10.8  ±  0.2 

-0.7 

-0.1 

0.0 

2MASS  J0501— 0010 

L4y 

14.3  ±  0.2 

13.2  ±0.2 

12.3  ±  0.2 

-1.0 

-0.7 

-0.7 

GL564c 

L4C 

12.5  ±  0.5 

12.0  ±  0.5 

11.0  ±0.5 

0.8 

0.4 

0.7 

2MASS  J0103+1935 

L6/3 

14.5  ±  0.4 

13.3  ±0.4 

12.5  ±  0.4 

-0.6 

-0.3 

-0.2 

HD203030B 

L7.5C 

15.0  ±0.6 

13.8  ±0.2 

13.1  ±0.2 

-0.8 

-0.5 

-0.5 

HN  Peg  B 

T2.5C 

14.5  ±0.1 

14.1  ±0.1 

13.8  ±0.1 

-0.1 

-0.1 

0.2 

Notes. 

a  A  values  are  calculated  from  the  MKO  Mjhk  values  from  the  polynomials  presented  in  Table  7.  Negative  values  indicate  underluminous  objects. 
b  2MASS  J0608— 2753  is  a  member  of  the  Pictoris  moving  group  as  described  in  Rice  et  al.  (2010). 
c  Objects  that  are  companions  to  young  (<^1  Gyr)  main-sequence  stars. 
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Figure  9.J  —  K  vs.  Mk  diagrams  for  L  and  T  dwarfs  with  the  evolutionary  models  of  Burrows  et  al.  (2006)  overplotted.  The  top  left  panel  shows  the  cloudless  model 
with  all  metallicities  and  the  remaining  three  panels  show  the  cloudy  model.  Varying  the  metallicity  from  subsolar  (0.0)  to  super-solar  (3.0)  and  gravity  from  low 
(log(g)  =  4.5)  to  high  (log(g)  =  5.5)  fits  the  L  and  T  dwarf  sequences  with  varying  degrees  of  accuracy. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


The  two  latest-type  L  dwarfs  in  our  sample  are  significantly 
redward  of  any  of  the  Burrows  et  al.  (2006)  predictions  but 
within  the  gravities  explored  on  the  cloudy  model  of  Saumon  & 
Marley  (2008).  Recent  work  by  Barman  et  al.  (201 1)  investigat¬ 
ing  HR8799b’s  red  near-IR  color  and  relatively  smooth  near-IR 
spectrum  concluded  that  thick  photospheric  dust  cloud  opacity 
could  explain  the  planets  observed  luminosity  and  color.  The 
latest-type  L  dwarfs  in  our  sample  may  be  higher  mass  analogs 
to  HR  8799b  (Barman  et  al.  2011). 

7.  SUBDWARFS 

There  are  12  ultracool  subdwarfs  with  parallaxes  including 
8  late-type  M  and  4  L  subdwarfs.  Figure  13  shows  the  near-IR 
absolute  magnitude  versus  spectral  type  diagrams  for  normal 
mid-type  M  through  late-type  L  dwarfs  (excluding  binaries  and 
low  surface  gravity  dwarfs)  with  the  subdwarfs  overplotted.  We 
show  all  subdwarfs  including  objects  with  absolute  magnitude 
uncertainty  >0.5  mag.  As  noted  in  Burgasser  et  al.  (2008d), 
the  L  subdwarfs  are  overluminous  in  Mj  but  shift  to  normal 
or  slightly  underluminous  by  MK.  This  has  been  attributed  to 
reduced  condensate  opacity,  as  evidenced  by  strong  TiO,  Ca  i, 
and  Ti  i  features  and  enhanced  collision-induced  H2  opacity  at 
the  K  band  (e.g.,  Ackerman  &  Marley  2001;  Tsuji  et  al.  1996; 
Burgasser  et  al.  2003a,  2007).  The  effect  is  not  as  pronounced  for 
the  late-type  M  subdwarfs  which,  with  the  exception  of  SSSPM 
J1 256— 1408,  appear  at  most  slightly  overluminous  in  Mj  and 
normal  or  underluminous  in  MH  and  MK. 


We  compare  photometry  for  the  two  L  subdwarfs 
2MASS  J0532+8246  and  2MASS  J1626+3925  in  Figures  8 
and  9.  In  Section  5.1,  we  find  that  early  and  mid- type  L  dwarfs 
are  best  fit  by  the  cloudy  tracks  using  both  Saumon  &  Marley 
(2008)  and  Burrows  et  al.  (2006)  evolutionary  models.  However, 
variations  in  the  metallicity  of  the  cloudy  tracks  in  Figure  9  do 
not  reproduce  the  colors  of  the  L  subdwarfs.  Instead,  one  must 
use  the  cloudless  tracks  supporting  the  idea  that  these  objects 
are  blue  and  overluminous  at  J  due  to  reduced  cloud  opacity. 

8.  KINEMATICS 

Combining  the  absolute  parallax  with  the  relative  proper 
motion  gives  the  tangential  velocity  (Vtan)  of  a  source  (see 
Table  2  for  Vtan  values  of  objects  studied  in  this  work).  As  our 
full  astrometric  sample  is  composed  of  objects  in  the  immediate 
solar  vicinity,  Vtan  values  and  their  dispersions  can  be  used  as  a 
rough  indicator  of  age.  In  general,  older  objects  will  have  had 
enough  time  to  interact  with  objects  in  the  Galactic  disk  and 
have  their  orbits  perturbed  while  younger  objects  will  retain  a 
motion  consistent  with  that  of  the  Galactic  disk  (i.e.,  co-moving 
with  their  nascent  cloud).  The  dispersion  of  a  population  is  more 
informative  than  individual  values  for  determining  ages. 

The  median  Vtan  and  crVdn  values  for  the  7 1  normal  and  unre¬ 
solved  binary  L  dwarfs  with  absolute  magnitude  uncertainties 
<0.5  are  27  km  s-1  and  20  km  s-1,  respectively.  For  the  50  nor¬ 
mal  and  unresolved  binary  T  dwarfs  we  find  similar  values  of 
31  km  s-1  and  20  km  s-1,  respectively.  Our  results  are  in  good 
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Figure  10.  7  —  K  vs.  Mk  diagram  for  L  and  T  dwarfs  with  the  evolutionary 
models  of  Saumon  &  Marley  (2008)  overplotted  to  demonstrate  the  best 
fit  for  L/T  transition  objects.  For  late-type  L  dwarfs  {Mk  <  13.0)  we  have 
overplotted  the  /se d  =  2  tracks  and  for  late-type  T  dwarfs  {Mk  >15.0)  we  have 
overplotted  the  /se d  =  4  tracks.  For  the  L/T  transition  we  created  a  hybrid  model 
between  the  two  by  adding  the  predicted  model  magnitudes  of  the  latter  to  the 
former  in  10%  increments  from  13.0  <  Mk  <  15.0. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


agreement  with  earlier  population  analyses  (Faherty  et  al.  2009; 
Schmidt  et  al.  2010).  Compared  to  the  kinematic  results  of  Vrba 
et  al.  (2004),  we  do  not  find  a  significant  difference  between 
the  kinematics  of  L  and  T  dwarfs.  Based  on  the  difference  in 
dispersion  between  the  L  and  T  dwarfs,  Vrba  et  al.  (2004)  con¬ 
cluded  from  their  much  smaller  sample  of  UCDs  that  the  L 
dwarfs  were  a  kinematically  younger  population  than  the  T 
dwarfs.  In  part,  their  conclusion  was  drawn  from  the  fact  that 
there  were  no  T  dwarfs  with  Etan  values  <  20  km  s  1 .  In  our 
larger  sample,  we  find  14  T  dwarfs  with  Vtan  values  <  20  km  s_1 . 


Spectral  Type 

Figure  11.  Spectral  type  vs.  absolute  magnitude  in  the  MKO  JHK  filters  for  late- 
type  M  through  mid-L  dwarfs.  Unfilled  circles  are  normal  dwarfs  with  parallax 
measurements.  Red  five  point  stars  and  filled  circles  are  intermediate  0 3 )  and 
low  (y)  surface  gravity  dwarfs,  and  purple  filled  circles  are  young  (<<Cl  Gyr) 
companions  to  nearby  stars.  In  Table  10,  we  report  the  difference  in  magnitude 
for  each  source  from  the  Mjhk  value  calculated  from  the  polynomials  in  Table  7. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


A  two-sided  K-S  test  on  the  L  and  T  velocities  yields  a  signifi¬ 
cant  probability  (p  ~  0.33)  that  the  L  and  T  dwarfs  in  our  sample 
have  identical  kinematics  hence  ages. 

We  isolated  the  low  surface  gravity  dwarfs  and  the  subdwarfs 
and  compared  their  kinematics  to  the  overall  sample  (note 
that  we  have  included  late-type  M  dwarfs  in  each  subset  but 
only  discuss  objects  whose  absolute  magnitude  uncertainties 
are  <0.5  mag).  The  former  have  significantly  smaller  Vtan 
values  and  tighter  dispersions  than  the  overall  population  and 
the  latter  significantly  larger  values.  The  median  Vtan  and  atan 
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Figure  12.  J  —  K  vs.  Mk  diagram  with  the  evolutionary  models  of  Saumon  & 
Marley  (2008;  top  panel)  and  Burrows  et  al  (2006;  bottom  panel)  overplotted 
along  with  four  low  surface  gravity  L  dwarfs.  The  log(g)  =  [4. 5, 5.0, 5. 5]  and 
/Sed  =  1,2  parameters  of  the  former  are  shown  as  are  the  log(g)  =  [4. 5, 5.0, 5.5], 
Metallicity=[0.0,3.0]  parameters  of  the  later. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

values  for  the  10  low  surface  gravity  dwarfs  are  10  km  s_l  and 
14  km  s_1 ,  respectively.  For  the  nine  subdwarfs,  the  median  Vtan 
and  crtan  values  are  241  km  s-1  and  68  km  s_1,  respectively.  The 
considerable  difference  in  values  for  each  subset  compared  to 
the  overall  population  further  confirms  expectations  that  they 
are  younger  (low  surface  gravity  objects)  and  older  (subdwarfs) 
than  the  overall  UCD  population. 

9.  BINARIES 

When  found  as  companions,  brown  dwarfs  are  primarily 
tightly  bound  (separations  <20  AU)  with  a  mass  ratio  close 
to  1  (e.g.,  Burgasser  et  al.  2003c;  Close  et  al.  2003).  As  such, 
near-equal  mass  brown  dwarf  binaries  can  be  identified  on  an 
H-R  diagram  by  their  overluminosity  compared  to  equivalent 
spectral  type  objects.  Moreover,  if  component  photometry  has 
been  acquired  one  can  investigate  the  properties  of  co-evolving 
systems  and  measure  dynamical  masses  (e.g.,  Dupuy  et  al.  2009, 
2008;  Konopacky  et  al.  2010).  Among  the  parallax  sample  listed 
in  Table  5  there  are  25  binaries  unresolved  in  2MASS.  These 
are  shown  in  Figure  3  as  filled  red  circles. 

There  were  four  sources  with  parallaxes  reported  in  this  work 
with  published  high  resolution  imaging  and  predicted  compo¬ 
nent  spectral  types  from  photometric  data  and/or  spectral  tem¬ 
plate  fitting:  2M0518— 2828,  2M1209-1004,  2M1404-3159, 
and  2M2052— 1609.  These  four  sources  were  used  in  the 
analysis  of  Section  4.1  and  their  system  properties  are  listed 
in  Table  6.  When  compared  to  the  absolute  magnitudes  calcu¬ 
lated  for  each  spectral  subtype  using  the  polynomial  reported  in 


Spectral  Type 

Figure  13.  Spectral  type  vs.  absolute  magnitude  in  the  MKO  JHK  filters  for 
late-type  M  through  L  dwarfs.  Unfilled  circles  are  normal  dwarfs  with  parallax 
measurements.  Blue  five  point  stars  are  subdwarfs  with  parallax  measurements. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Table  7,  both  components  of  2M 1209— 1004,  2M2052— 1609, 
and  2M0518— 2828  fit  within  uncertainties  on  the  brown  dwarf 
main  sequence.  While  the  primary  for  2M1404— 3159  also  fits 
within  the  brown  dwarf  main  sequence,  the  T5  secondary  is 
~1  mag  underluminous  in  MJHK. 

Burgasser  et  al.  (2010a)  used  a  template  fitting  technique 
to  identify  17  new  L-T  transition  binaries  and  three  of  the 
suspected  binaries  have  parallaxes  reported  in  this  work: 
SDSS151 1+0607,  2MASSJ0949— 1545,  and  SDSS 1207+0244. 
The  first  was  regarded  as  a  strong  binary  candidate  while  the 
latter  two  are  weak  binary  candidates.  Examining  the  H-R  dia¬ 
gram  in  Figure  3  we  find  SDSS  15 11+0607  to  be  nearly  1  mag 
overluminous  in  MJHK.  C.  R.  Gelino  et  al.  (in  preparation) 
have  followed  up  with  Adaptive  optics  and  resolved  the  two 
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components.  We  show  both  components  in  Figure  5  and  list  the 
individual  magnitudes  in  Table  6.  The  proposed  L5.5  and  T5 
components  fit  well  within  the  brown  dwarf  main  sequence. 

Using  new  parallax  measurements,  2MASSJ0949— 1545  and 
SDSS 1207+0244  fit  well  on  the  H-R  diagram  as  single  objects. 
If  decomposed  into  the  binary  components  proposed  by  Bur- 
gasser  et  al.  (2010a),  the  primary  and  secondary  sources  would 
be  >  1  mag  underluminous  in  multiple  bands.  We  conclude  that 
these  sources  are  best  treated  as  single  sources. 

10.  CONCLUSIONS 

We  have  measured  parallaxes  for  11  M,  32  L,  and  27  T 
dwarfs  in  the  local  solar  neighborhood.  Nine  calibrator  stars 
were  included  in  the  sample  to  verify  the  reliability  of  our 
pipeline.  The  focus  of  this  project  was  on  low  surface  gravity 
dwarfs,  L/T  transition  objects,  and  late-type  T  dwarfs  within 
20  pc  of  the  Sun.  The  70  new  parallaxes  significantly  increase  the 
number  of  brown  dwarfs  with  accurate  distance  measurements. 

We  combined  our  sample  with  115  literature  measurements 
and  used  the  full  sample  to  re-define  color-magnitude  and 
spectrophotometric  diagrams  in  JHK.  Adding  decomposed 
L/T  transition  binaries  we  find  no  reason  to  split  spectral-type/ 
absolute  magnitude  polynomials  into  “bright”  and  “faint”  trends 
to  account  for  unresolved  binarity  as  has  been  done  in  the  past. 
Isolating  T0-T4  dwarfs  to  investigate  the  extent  of  the  L/T 
transition  brightening  we  find  that  there  is  a  [1.2-1. 4]  mag 
difference  at  /,  [0.3-0. 5]  mag  difference  at  H,  and  a  plateau 
or  dimming  of  [—0.2  to  —0.3]  mag  at  K.  In  agreement  with 
flux  reversal  binary  studies,  this  confirms  the  brightening — and 
the  physical  mechanism  that  drives  it — as  an  intrinsic  feature  of 
brown  dwarf  atmospheric  evolution. 

We  compared  the  /  —  K  versus  MK  data  for  the  full  parallax 
sample  to  the  evolutionary  models  of  Saumon  &  Marley  (2008) 
and  Burrows  et  al.  (2006).  The  /se d  =  1,2  parameters  best  fit  the 
L  dwarf  sequence  and  /se d  =  4  (corresponding  to  a  very  thin 
cloud  layer)  best  fit  the  late-type  T  dwarf  sequence  using  the 
Saumon  &  Marley  (2008)  models.  The  cloud  model  with  varying 
gravity  and  metallicity  reproduces  the  L  dwarfs  and  the  clear 
model  with  similar  variations  fits  the  T  dwarfs  using  the  Burrows 
et  al.  (2006)  models.  However,  comparisons  of  both  models  to 
empirical  data  show  significant  red  or  potentially  “ultra-cloudy” 
L  dwarf  outliers.  Similarly  there  is  significant  scatter  seen  in  the 
latest  type  T  dwarfs  that  is  unaccounted  for  in  the  clear  and 
/sed  =  4  models  indicating  that  condensate  clouds  continue 
to  play  a  role  in  the  photospheres  of  some  low-temperature 
brown  dwarfs.  Investigations  of  individual  objects  in  the  “ultra- 
cloudy”  sample  reveals  objects  which  are  young/moderately 
young,  implying  a  correlation  between  youth  and  enhanced 
photospheric  dust. 

No  single  /sed  parameter,  gravity,  nor  metallicity  track  in 
the  evolutionary  models  can  account  for  the  L/T  transition 
objects.  However,  a  hybrid  model  that  smoothly  transitions  from 
/sed  =  2  (the  best  fit  for  late-type  L  dwarfs)  to  /se d  =  4  (the  best 
fit  for  late-type  T  dwarfs)  at  a  near-constant  reff  =  1200  zb  100  K 
encompasses  the  majority  of  T0-T4  dwarfs.  This  temperature 
range  is  consistent  with  recent  toy  and  sophisticated  hybrid 
models  but  demonstrates  that  the  range  of  spectral  subtypes  for 
which  the  temperature  plateau  is  applicable  is  narrower  than 
previously  suspected  (T0-T4  as  opposed  to  L7-T4). 

The  low  surface  gravity  objects  with  parallax  measurements 
in  this  work  are  not  explained  by  varying  gravity  in  the 
evolutionary  models.  Rather,  they  appear  to  be  underluminous 
compared  to  model  color-magnitude  diagrams  and  for  their  low 


gravity  spectral  type.  Among  the  10  low  surface  gravity  or  young 
companion  L  dwarfs  investigated,  80%  appear  [0.2-1. 0]  mag 
underluminous  for  their  spectral  type  in  /,  H,  and/or  K.  Possible 
explanations  for  their  underluminosity  are  that  (1)  the  low- 
gravity  and  field  dwarf  spectral  classification  schemes  are  on 
different  temperature  scales  with  the  low  surface  gravity  objects 
intrinsically  cooler  than  field  age  objects  of  the  same  type; 
and/or  (2)  young  objects  could  have  dustier  photospheres  than 
field-aged  objects  making  them  appear  both  fainter  and  redder. 

A  kinematic  analysis  of  the  astrometric  sample  reveals  similar 
velocity  dispersions  between  the  L  and  T  dwarf  populations.  A 
two-sided  K-S  test  verifies  that  the  two  kinematic  distributions 
are  likely  drawn  from  the  same  population  and  hence  have 
similar  age  distributions.  The  low  surface  gravity  and  subdwarf 
samples  have  distinctly  different  velocity  dispersions  and  are 
likely  significantly  younger  and  older  (respectively)  than  normal 
objects. 
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